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1.  Introduction  and  Summary 


This  report  documents  work  performed  for  the  Air  Force  Geophysics 
Laboratory  (AFGL)  and  the  Air  Force  Space  Division  (AFSD)  on  the  subject 
of  optical  emissions  from  the  ionosphere  and  their  relationship  to 
electron  density  profiles  (EDPs).  The  work  spans  the  period  from 
February  1983  to  September  1984,  with  actual  completion  several  months 
prior  to  the  termination  date.  Considerable  documentation  was  produced 
under  this  contract; 

•  Quarterly  reports  (6/83,  9/83,  12/83,  and  3/84) 

•  Technical  report  to  AFSD  (Strickland  et  al.,  1984a) 

•  This  final  report 

•  Publications 

-  Strickland  et  al.  (1984b) 

-  Daniell  et  al .  (1984) 

-  Daniell  and  Strickland  (1985). 

Copies  of  the  publications  appear  as  appendices  to  this  report. 


Prior  contract  work  concentrated  on  the  study  of  auroral  processes 
using  first-principles  models  and  applying  these  to  making  predictions 
and  analyzing  rocket  and  satellite  data.  Interest  focussed  on  the 
continuous  aurora  to  better  understand  chemical  and  radiation  processes 
and  determine  how  useful  satellite  optical  data  can  be  in  deducing  the 
auroral  E  EDP  (see  Strickland  et  al.,  1983a).  Work  during  the  present 
contract  expanded  to  include  investigation  of  ionospheric  processes  in 
the  mid-latitude  nighttime  and  low  to  mid-latitude  daytime  ionospheres. 
Much  of  our  time  was  directed  to  the  nighttime  problem  which  is  briefly 
outlined  in  the  next  section,  and  has  been  documented  in  the  above 
referenced  AFSD  report  and  in  the  paper  by  Daniell  et  al .  (1984).  We 
have  determined  in  quantitative  terms  how  well  one  can  expect  to  deduce 
the  nighttime  EDP  under  quiescent  plasma  conditions  from  the  weak  atomic 
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oxygen  emissions  resulting  from  recombination  of  0^^  and  0^  with  ambient 
electrons.  Present  state-of-the-art  detectors  recording  data  from  a 
satellite  such  as  DMbP  will  significantly  improve  the  Air  Force's  present 
near  real-time  specification  of  this  EOP. 

We  may  summarize  our  progress  and  findings  as  follows: 

1.  Two  nighttime  codes  have  been  developed  and  applied.  One 
specifies  the  nighttime  EDP  given  nadir  intensities  of  01  6300 
and  01  1356  A.  The  other  does  the  reverse  calculation.  Both 
codes  use  an  empirical  (Chapman)  model  to  characterize  the  shape 
of  the  hUP. 

2.  A  literature  survey  of  emissions  and  the  EDP  under  the  nighttime 
conditions  was  conducted.  The  survey,  plus  application  of  the 
above  codes,  have  provided  good  morphological  information  on 
these  quantities,  giving  their  behavior  as  functions  of  local 
time,  season,  solar  activity,  and  latitude. 

3.  The  nighttime  background  to  be  found  in  the  bandpass  of  a  6300  A 
photometer  was  investigated.  Its  contribution  to  the  6300  A 
signal  was  determined  for  various  tilting  filter  configurations 
and  over  the  full  range  of  lunar  luminosities.  It  was 
determined  that  a  2A  tilting  filter  photometer  should  be  able  to 
effectively  determine  the  true  01  6300  A  intensity  within 
integration  time  constraints. 

4.  An  investigation  was  carried  out  which  provided  instrument 
specifications  needed  to  adaquately  measure  optical  intensities 
at  night,  in  the  day,  and  for  the  continuous  aurora.  The 
specifications  for  day  and  auroral  observations  applied  to  a 
spectrometer  while  those  for  nighttime  observations  applied  to 
photometers.  Section  4  reports  our  findings. 

5.  An  investigation  of  auroral  EUV  emissions  was  conducted  to 
determine  their  variability  and  possible  use  in  deducing  EDPs 
(Strickland  et  al.,  1983b).  This  was  documented  in  the  second 
quarterly  report. 
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6.  A  code  was  developed  which  determines  limb  intensity  profiles 
for  optically  thick  or  thin  emissions  with  or  without  pure 
absorption.  The  algorithim  is  efficient  and  guarantees  a  high 
degree  of  accuracy  since  the  integrand  is  investigated  prior  to 
integration.  This  is  documented  in  the  first  quarterly  report. 

7.  Modeling  of  the  daytime  EDP  and  selected  optical  emissions  was 
conducted  to  compare  with  rocket  and  ground-based  sounder  data 
and  with  independent  modeling  results  being  obtained  at  AFGL  by 
Decker  and  Jasperse.  See  Section  5  for  a  discussion  of  the 
daytime  problem. 

8.  Further  analysis  was  carried  out  on  the  optical  and  EDP  data 
from  the  AFGL  auroral  E  rocket  program.  Good  agreement  was 
achieved  with  the  available  data.  Details  are  available  in 
Section  5  of  the  AFSD  report.  Publication  of  the  analysis  is 
planned  in  1985. 

9.  An  ongoing  search  was  performed  for  experiments  (rocket, 
satellite,  and  ground  based)  yielding  simultaneous  EDP  and 
optical  data  in  the  three  global  regions  cited  above.  The  most 
notable  example  found  is  discussed  in  Section  5. 

In  summary,  several  investigations  were  conducted  with  emphasis  on 
determining  how  well  optical  data  can  be  used  to  deduce  the  EDP  in  the 
nighttime,  daytime,  and  auroral  ionospheres.  The  probleii'  was  limited  to 
the  E  region  in  the  case  of  the  aurora.  Both  the  E  and  F  regions  were  of 
interest  for  the  other  two  regions.  We  have  much  better  quantitative 
information  on  the  use  of  this  concept  than  when  this  work  began.  It 
looks  particularly  attractive  at  night  where  the  emissions  provide  a 
direct  signature  of  the  EDP.  Much  more  work  is  still  needed,  especially 
with  regard  to  F  region  global  behavior. 

A  satellite  experiment  is  also  needed  to  definitively  test  our  models 
being  used  to  relate  optical  intensities  to  EDPs.  Such  an  experiment 
will  probably  begin  in  1989  under  the  direction  of  Meier  and  colleagues 
at  the  Naval  Research  Laboratory.  The  instrument  is  called  RAIDS  (Remote 
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Atmospheric  and  Ionospheric  Detection  System).  It  will  measure 
intensities  from  the  EUV  to  the  visible  ror  various  viewing  geometries 
(including  nadir  and  limb  viewing).  Simultaneous  EDP  measurements  will 
occur  when  the  satellite  is  passing  over  ground-based  ionosonde  stations. 
Huffman  and  colleagues  at  AF6L  are  also  beginning  work  on  an  operational 
system  to  fly  on  DMSP  in  the  late  1980s  and  early  1990s.  This  is  the 
system  for  which  we  have  been  conducting  feasibility  studies  over  the 
pase  several  years.  Nadir  observed  intensities  and  in  situ  measurements 
of  the  electron  density  and  temperature  are  being  planned  with 
simultaneous  measurements  of  the  EDP  by  ground-based  ionosondes  whenever 
DMSP  is  approximately  overhead.  The  RAIDS  experiment,  which  is  research 
oriented,  should  provide  valuable  information  to  us  prior  to  the  start  of 
Huffman's  operational  system. 
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2.  Overview 

We  have  included  this  section  to  address  the  following  items: 

•  Larger  Air  Force  interests 

•  Categories  under  which  the  completed  work  may  be  broken  down 

•  Documentation  in  this  report  versus  total  documentation  under 
this  contract 

This  work  is  motivated  by  the  need  to  improve  the  present  near  real 
time  specification  of  the  EDP  on  a  global  scale  from  the  lower  E  region 
to  the  upper  F  region.  Systems  in  need  of  this  information  include  HF 
communications  systems,  over-the-horizon  (OTH)  radar  systems,  SPACETRACK 
radar  systems,  and  classified  systems.  Satellites  are  the  measuring 
platforms  needed  for  global  coverage  and  offer  two  options  for  monitoring 
the  EDP:  active  or  remote  sensing.  Active  sensing  requires  the  use  of  an 
ionosonde  which  yields  the  EDP  maximum.  Passive  sensing  requires  the  use 
of  photometers  and/or  spectrometers  capable  of  measurinr  radiations  over 
a  wide  intensity  range.  These  intensities  can  be  related  to  EDP 
characteristics  as  well  as  the  dominant  neutral  densities.  The  key  issue 
here  relates  to  the  size  of  error  bars  placed  on  the  deduced  EDP  and 
whether  it  is  small  enough  to  significantly  improve  present  global  EDP 
specification. 

We  have  been  investigating  the  passive  sensing  problem  for  the  past 
few  years.  Prior  to  this  contract,  work  concentrated  primarily  on  the 
aurora  -  in  particular,  the  continuous  aurora  known  for  its  stability 
over  tens  of  minutes  (see,  e.g.,  Whalen,  1981).  Good  examples  of  studies 
undertaken  may  be  seen  in  the  papers  by  Strickland  et  al .  (1983a)  and 
Daniell  and  Strickland  (1985).  We  have  examined  in  detail  how  a  variety 
of  emissions  from  the  EUV  to  the  visible  along  with  the  EDP  change  as 
functions  of  the  precipitating  electron  spectrum  and  atmospheric 
conditions.  We  found  that  selected  emission  efficiencies  changed 
noticeably,  showing  potential  for  monitoring  electron  precipitation 
character  1  St  ics  which  in  turn  can  yield  the  EDP  through  modeling.  This 
work  concentrated  on  FOP  specification  for  the  F  region. 


In  this  contract,  emphasis  has  shifted  from  the  auroral  ionosphere  to 
the  mid-latitude  nighttime  and  daytime  ionospheres.  Part  of  this  was 
caused  by  the  discovery  of  strong  clutter  in  OTH  returns  from  the  auroral 
ionosphere  during  radar  experiments  conducted  by  the  Air  Force.  The 
clutter  is  produced  by  plasma  structure  arising  from  instabilities  and 
reduces  the  potential  of  OTH  radar  as  an  effective  tool  for  detecting 
aircraft  and  missiles  at  high  northern  latitudes.  Beyond  this,  there  was 
a  need  to  begin  examining  the  optical  concept  elsewhere  since  interest  is 
on  a  global  scale.  As  noted  in  the  introduction,  most  of  our  effort  has 
been  directed  to  the  nighttime  ionosphere.  Some  work  has  continued  on 
auroral  problems,  particularly  on  analysis  of  data  from  the  AFGL  auroral 
E  rocket  program.  Some  work  has  also  been  carried  outon  the  daytime 
problem,  mainly  in  collaboration  with  personnel. 

Table  1  shows  a  breakdown  into  eight  categories  of  the  work 
undertaken  in  this  contract.  We  see  that  most  catagories  apply  to  the 
nighttime  investigation.  A  brief  description  of  this  investigation 
follows.  We  began  by  developing  two  codes  which  relate  the  01  6300  A  and 

01  1356  A  emissions  (there  are  other  01  lines  such  as  1304  A  which  could 

serve  the  same  purpose  as  1356  A)  to  the  EDP.  One  code  calculated  an  EOP 

(assuming  a  Chapman  dependence)  from  the  emissions  while  the  other  does 

the  reverse  calculation  given  a  particular  EDP.  These  codes  were  used  to 
predict  the  morphological  behavior  of  these  emissions  and  to  examine  the 
sensitivity  of  the  EDP  to  variaions  in  either  emission.  Another 
important  aspect  of  the  work  was  to  characterize  the  background  around  01 
5300  A  and  do  statistical  studies  to  determine  how  long  this  feature  must 
be  observed  to  extract  a  good  6300  A  signal.  The  assumed  detection 
system  was  a  tilting  filter  photometer  with  various  bandpasses.  Such  a 
photometer  is  capable  of  scanning  in  wavelength,  which  enables  the  true 
signal  to  be  separated  from  the  background  signal.  Detectability  was 
investigated  for  various  levels  of  background  and  01  6300  A  intensity. 

The  most  detailed  documentation  of  the  nighttime  investigation  appears  in 
our  recent  AFSD  report  (Strickland  et  al.,  1984a). 

We  complete  this  discussion  with  the  third  topic  listed  above.  Most 
of  the  work  reported  in  the  following  sections  has  already  been 
documented  elsewhere.  The  next  section  provides  a  synopsis  of  the  major 


f 


p. 

k  • 


n  • 

p- 


r-;- 


vnv  ■  ^  u  ■  u  ■  u  •  L' » WWW  V.^"  **  v*  T" 


document  to  be  generated  in  this  contract  (the  report  just  referenced 
above).  The  appendices  contain  papers  either  published  or  under  review. 
Section  4  contains  a  revised  internal  report  originally  prepared  for 
Huffman  at  AFGL.  The  subject  deals  with  instrumentation  and  gives 
requirements  for  the  detection  of  key  optical  features  in  the  nighttime, 
daytime,  and  auroral  ionospheres.  This  represents  a  beginning  to  an 
obviously  critical  part  of  the  program,  if  the  concept  is  to  lead  to  an 
operational  system.  Section  5  discusses  a  small  effort  directed  to 
analysing  data  from  a  daytime  rocket  experiment.  The  experiment  provided 
the  opportunity  to  analyze  coincident  optical  and  EDP  data.  Such 
analyses  are  critical  to  placing  error  bars  on  EDPs  to  be  deduced  from 
optical  data. 

Table  1.  Breakdown  of  contract  work  by  region  of  the  ionosphere 


Categories: 


NIGHT  DAY  AURORA 


Code  development 
First  principles  modeling 
Empirical  Modeling 
Literature  Survey 
Gathering  data 
Analysing  Data 
Detectabi 1 i ty 
Backgrounds 

Instrumental  characteristics 
Statistics 
Reporting 


X 

X 

X 

X 

X  X 

X  X 

X 

X  X 

X 

X  X 


X 


X 


X 
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3.  The  Proposed  UV  Ionospheric  Monitoring  System 


Much  of  the  work  we  performed  under  this  contract  was  described  in 
Strickland  et  al,  (1984a).  Rather  than  repeat  that  material  here,  we 
have  chosen  to  provide  a  synopsis  based  on  the  executive  summary  of  that 
report.  The  effort  was  directed  toward  developing  a  system  concept  for  a 
satellite-borne  optical  sensor  and  the  associated  data  analysis  software 
that  would  be  capable  of  providing  electron  density  profiles  over  much  of 
the  globe  and  in  near  real  time. 


3.1  The  System  Concept 

The  system  concept  includes  (a)  developing  and  installing  a  UV  sensor 
on  the  OMSP  satellite  and  developing  an  associated  automatic  data 
processing  (software)  system,  (b)  using  the  optical  data  with  DMSP 
electron  density  and  temperature  data  to  deduce  the  near  real-time  EDP 
near  the  satellite  orbital  plane,  and  (c)  transmitting  the  EDP  data  for 
use  with  ground-based  ionosonde  data  and  GPS  total  electron  content  (TEC) 
data  to  specify  the  global  EDP  for  system  users.  The  ionosonde  and  TEC 
data  would  serve  two  purposes:  to  augment  the  EDP  data  base  generated  by 
the  UV  data,  and  to  refine  the  UV-based  EDP  where  coincident  data  exist. 


3.2  The  Ionospheric  Subregions 

We  have  considered  the  following  ionospheric  subregions;  (a)  the 
daytime  low-  to  mid-latitude  ionosphere  from  90  to  1000  km,  (b)  the 
nighttime  mid-latitude  ionosphere  from  250  to  1000  km,  and  (c)  the 
auroral  E  layer  from  90  to  200  km  for  undisturbed  conditions.  The 

spatial  resolution  considered  for  the  daytime  and  nighttime  EDP  is  at 

least  one  vertical  profile  for  each  square  500  km  on  a  side  and  for  the 
auroral  E  layer,  at  least  one  vertical  profile  for  each  50  km  on  a  side. 

Because  of  the  high  frequency  of  occurance  of  irregularities  or  of  highly 

variable  transport  conditions,  the  polar  caps,  the  cusp,  the  auroral  F 
region,  and  the  equatorial  nighttime  regiori  were  excluded  from 
consideration. 
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3.3  Use  of  UV  Measurements 


(a)  In  the  daytime  low-  to  mid-latitude  ionosphere  there  are  no 
useful  optical  emission  features  which,  when  measured  from  the  DMSh 
altitude  (  840  km),  give  a  signal  directly  dependent  upon  the  electron 
density  profile.  Therefore,  we  propose  an  indirect  method  in  which  we 
measure  the  emission  feature  at  1356  A  and  at  least  one  of  the 
'.yman-Birge-Hopf ield  (LBH)  bands  to  deduce  the  0  to  density  ratio  and 
the  absolute  scale  of  the  solar  flux.  This  information  and  a  neutral 
wind  model  are  used  to  calculate  the  EDP  from  first  principles  using  time 
dependent  ion  continuity  equations.  Basing  the  specified  EDP  on  the 
deduced  key  parameters  (0  to  density  ratio  and  scaling  factor  for  the 
solar  flux)  significantly  improves  the  accuracy  of  the  EDP  over  a  result 
based  on  empirical  values  of  these  parameters. 

(b)  In  the  nighttime  mid-latitude  ionosphere  there  are  at  least  two 
usable  emission  features  which  give  a  signal  that  is  directly  dependent 
on  the  EDP.  The  emission  features  are  the  atomic  oxygen  lines  at  1356  A 
and  6300  A.  These  features  are  most  sensitive  to  the  electron  content 
near  and  below  the  peak  of  the  EDP.  To  extract  the  mos'  accurate  EDP 
from  these  observations,  we  are  investigating  two  approaches:  (i)  a  time 
dependent  transport  approach  in  which  the  neutral  winds  and  possibly  the 
electric  fields  are  deduced  from  the  optical  emission  measurements,  and 
(ii)  use  of  an  empirical  model  in  which  the  shape  of  the  F  layer  is  fixed 
and  its  location  in  altitude  and  its  absolute  value  are  determined  by  the 
emission  features.  This  would  produce  F  layer  profiles  from  about  ?50  to 
1000  km. 

(c)  For  the  E  layer  of  the  diffuse  aurora  we  can  use  the  emission 
feature  at  1356  A  and  one  or  more  of  the  LBH  bands  to  determine  the 
hardness  and  energy  flux  of  the  incident  auroral  electrons  which  produce 
the  ionization.  The  EDP  is  then  determined  from  90  to  about  200  km  using 
this  incident  electron  information  in  our  electron  tr ansport/chemi stry 
model . 
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3.4  Uncertainties 


We  have  found  that  the  proposed  observables  can  provide  a  large 
fraction  of  the  required  90-1000  km  global  coverage  (see  Table  2),  as 
well  as  the  bonus  of  sensing  a  broad  area  (rather  than  a  single  line) 
under  the  sub-orbital  track.  The  feasibility  of  this  approach  rests  upon 
our  own  recent  advances  in  state-of-the-art  capabilities  simultaneously 
pressed  into  parallel  areas:  hardware  development  of  UV  sensors  capable 
of  measuring  complete  global  day  and  nighttime  emissions,  and  physical 
modeling  of  ionospheric  processes  strongly  coupled  to  selected  optical 
emissions.  The  hardware  development  has  demonstrated  success  on  both  a 
shuttle  flight  and  a  recent  satellite  project;  the  software  development 
has  been  successfully  tested  against  field  experimental  data. 

AVAILABLE  OPTIONS 


SUBREGIONS 

UV 

TOPSIDE 

SOUNDER 

X-RAY 

$  Daytime  low-  to  mid-latitude 

Bottoms ide 

Yes 

No 

No 

Topside 

Yes 

Yes 

No 

•  Nighttime  low-  to  mid-latitude 

Bottoms ide 

Yes 

No 

No 

Topside 

Yes 

Yes 

No 

•  Auroral  zone 

E  Layer 

Yes 

No 

Yes 

D  Layer 

No 

No 

Yes 

Table  2.  Capabilities  of  Three  Methods  for  Remote  Sensing  of  the 
Ionosphere 

Very  few  data  exist  for  quantitative  determination  of  error  bars  for 
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this  EDP  system  concept  due  to  scarcity  of  coincident  optical  and  EDP 
measurements.  For  the  present,  the  expected  EDP  improvement  must  be 
discussed  in  relative  terms  rather  than  in  absolute  terms.  Nevertheless, 
meaningful  statements  can  be  made. 

Errors  in  our  ab-initio  daytime  model  are  comparable  to  those  of  the 
other  state-of-the-art  models  for  known  conditions  (order  +  30%).  The 
factors  of  prime  importance,  as  noted  in  Section  1.3,  are  the  solar  EUV 
flux,  the  O/Np  density  ratio,  and  the  neutral  wind.  Errors  in  excess  of 
a  factor  of  2  can  be  expected  in  the  calculated  EDP  without  direct 
knowledge  of  these  factors.  Optical  observations  will  provide  the  needed 
knowledge  for  the  first  two  factors  and  should  lead  to  acceptable 
accuracy  in  the  E  and  lower  F  regions.  Near  and  above  the  peak  of  the 
EDP,  transport  effects  are  important.  Further,  overall  errors  should  be 
signif icantly  reduced  with  the  aid  of  ionosonde  and  TEC  data  cited 
previously.  Such  data  allow  spo',  checking  of  the  UV-deduced  EDP  through 
which  adjustments  can  be  made  f o  •  such  factors  as  optical  sensor 
degradation. 

Errors  in  deduced  nighttime  :i  profiles  arise  from  /ariable  transport 
and  uncertainties  in  temperature  and  light  ion  composition.  The  major 
effects  of  variable  transport  are  changes  in  the  overall  magnitude  of  the 
EDP  and  location  of  the  EDP  maximum.  A  less  important  effect  is  a  change 
in  the  EDP  shape.  The  proposed  observables  directly  specify  both  the 
electron  concentration  near  the  peak  and  the  height  of  the  peak,  thus 
eliminating  most  of  the  uncertainty  caused  by  variable  transport.  The 
less  important  shape  change  remains  a  source  of  EDP  error  for  this 
approach,  being  perhaps  20  to  30  percent  at  times  of  strong  transport. 
Topside  effects  due  to  temperature  and  light  ion  composition  are  likely 
to  be  significant.  The  optimum  approach  to  merging  the  information  from 
optical  emissions  with  information  from  in  situ  sensors  remains  to  be 
determined. 

In  the  auroral  E  layer,  for  arc-free  conditions,  determination  of  the 
EDP  from  90  to  200  km  should  be  achievable  to  within  +  40%  by  the  method 
proposed  in  this  report. 
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4.  Instrument  Requirements  for  a  Global  Ionospheric  Monitoring  System 

4.1  Introduction 

The  system  concept  described  in  the  preceding  section  involves  the 
use  of  optical  data,  in  situ  plasma  data,  and  a  software  data  analysis 
system.  We  have  made  a  considerable  effort  to  determine  the 
characteristics  of  the  optical  sensors  that  are  necessary  for  determining 
the  EDP  to  the  described  accuracy.  The  system  is  intended  to  monitor 
three  regions: 

(1)  daytime  low-  to  mid- latitude  ionosphere  from  90  to  1000  km 

(2)  nighttime  mid-latitude  ionosphere  from  250  to  1000  km,  and 

(3)  the  auroral  E  layer  from  90  to  200  km. 

Monitoring  is  to  be  done  under  quiescent  plasma  conditions.  The  spatial 
resolution  for  the  day  and  night  EDP  is  to  be  at  least  one  vertical 
profile  for  each  square  500  km  on  a  side  and,  for  the  auroral  E  layer,  at 
least  one  profile  for  each  square  50  km  on  a  side. 

The  system  is  intended  for  use  on  DMSP  satellites,  which  are  placed 
in  sun  synchronous  orbits,  that  is,  orbits  which  always  cross  the  equator 
at  the  same  local  time.  There  are  two  such  orbits  in  use.  One  is 
approximately  noon -midnight;  the  other  is  approximately  dawn-dusk.  The 
noon-midnight  orbit  is  shown  in  Figure  1.  It  crosses  the  equator  at 
local  times  of  10:30  am  and  10:30  pm.  The  dawn-dusk  orbit  is  also 
displaced  from  its  nominal  local  time  and  actually  crosses  the  equator 
before  dawn  and  before  dusk. 

4.2  Optical  Emissions 

Table  3  lists  some  of  the  candidate  optical  features  and  nominal 
column  emission  rates  for  the  three  global  regions  of  interest.  The 
absence  of  a  value  indicates  that  the  given  feature  is  not  useful  for  EDP 
monitoring  in  that  region.  The  nighttime  values  are  very  small  and  have 
important  implications  for  system  design.  However,  these  emissions  are 
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Typical  DMSP  Orbit 


i 

Noon 

Midnight 


Figure  1.  A  typical  DMSP  noon-midnight  orbit.  Note  that  the  orbit  actually  crosses  the 
_ equator  at  1030  and  2230  hr.  The  tick  marks  on  the  orbit  track  are  two  minutes 


very  attractive  since  they  provide  direct  signatures  of  the  EDP. 
Photometers  will  be  needed  to  accurately  measure  them.  Emissions 
in  the  day  and  during  auroral  activity  are  much  brighter.  A  UV 
spectrometer  will  be  capable  of  monitoring  them  on  time  scales  of 
interest.  The  capability  should  be  present  to  produce  images  (as  was 
done  on  HILAT  with  the  AIM  instrument). 

Monitoring  the  6300  A  emission  at  night  presents  special  challenges 
beyond  its  inherent  weakness.  Part  of  the  problem  can  be  seen  in  Figure 
2.  The  spectrum  comes  from  Broadfoot  and  Kendall  (1968)  and  shows  the 
presence  of  other  emissions.  In  particular,  unwanted  emission  comes  from 
OH  rotational  lines  and  from  backscattered  astronomical  light.  In 
addition,  a  nadir  viewing  instrument  will  register  photons  '"‘fleeted  from 
the  ground  or  clouds  as  well  as  photons  coming  directly  from  the  emitting 
layer.  We  propose  the  use  of  a  tilting  filter  photometer  or  equivalent 
system  to  monitor  both  the  signal  and  background. 

Examples  of  partial  FUV  spectra  for  day  and  auroral  conditions  are 
shown  in  Figure  3.  These  were  calculated  for  the  conditions  shown  and 
provided  intensities  for  additional  features  not  shown  in  Table  3.  We 
observe  variability  with  changing  conditions,  which  is  the  key  to  using 
the  proposed  features  for  EDP  monitoring. 


Feature 

Day 

Mid-latitude 

Night 

Auroral 

E  Layer 

1356 

100  -  3000 

.1-4 

100 

-  1000 

1304 

1000  -  20,000 

.2  -  10 

1000 

-  10,000 

911 

- 

.05  -  2 

- 

1273 

5-50 

- 

10 

-  50 

1325 

10  -  100 

- 

20 

-  100 

1383 

20  -  200 

- 

30 

-  150 

6300 

2  -  100 

• 

Table  3.  Nominal  Ranges  of  Column  Emission  Rates  in  Rayleighs  as 
observed  from  DMSP  Satellite  Altitude.  (Lack  of  entry  indicates  feature 
is  not  to  be  used  in  given  region.) 


nomical  background  (Broadfoot  and  Kendall,  1980). 
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4,3  The  Monitoring  System 


The  ideal  optical  monitoring  system  would  use  the  same  instrument  to 
monitor  all  three  regions  of  interest.  Unfortunately,  the  different 
characteristics  of  the  dayglow,  nightglow,  and  aurora  cause  difficulties 
in  designing  an  instrument  that  can  make  acceptable  measurements  in  all 
three  environments.  A  high  resolution  instrument  is  reguired  for 
measurements  of  the  important  UV  spectral  features  in  the  dayglow  and 
aurora.  These  features  include  01  1356  A  and  LBH  bands  at  1383  A, 

1672  A,  and  1768  A.  It  would  be  highly  desirable  to  monitor  the  1366  A 
feature  and  one  or  more  of  the  LBH  bands  simultaneously.  This  will  be 
done  by  the  AFGL  UV  imager  on  the  Polar  Bear  Satellite.  At  night, 
instrument  sensitivity  is  more  important  than  resolution.  The  UV 
features  of  interest  are  01  1356  A  and  01  1304  A,  which  may  be  observed 
together  by  a  wide  band  photometer.  In  addition,  the  01  6300  A  line  must 
also  be  monitored,  and  that  monitoring  requires  a  separate  photometer. 
Thus,  the  complete  system  should  consist  of  an  imaging  UV  spectrometer 
capable  of  observing  two  or  more  features  simultaneously,  a  UV 
photometer,  and  a  visible  photometer. 

To  quantify  the  measurement  requirements,  we  introduce  R,  the 
instrument  response  function  (or  sensitivity),  which  applies  to  any 
optical  system.  It  is 

R  =  (10®/4it)  A^tQT  counts  Rayleigh"^  sec"^ 


where 
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A  is  the  area  of  entrance  aperture,  (cm  ) 
is  the  viewing  solid  angle,  (sr) 

0  is  the  quantum  efficiency,  and 
T  is  the  transmission  function 

The  scaling  factor  of  10^/4  it  allows  R  to  be  expressed  in  terms  of 
Rayleighs. 


In  the  following  discussion,  we  specify  the  minimum  value  of  R  that 
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an  instrument  must  have  to  make  the  measurements  with  the  required 
accuracy.  We  do  not  specify  values  for  aperture,  solid  angle,  or  other 
parameters,  because  they  are  also  constrained  by  engineering  requirements 
that  we  have  not  addressed.  We  have  found  that  in  some  cases  the 
required  instrument  response  is  much  greater  than  the  response  that  can 
be  readily  obtained.  For  reference,  we  note  that  Huffman  et  al .  (1980) 
describe  a  UV  spectrometer  with  R  =  55  counts  R”^  s”^.  (This  instrument 
was  flown  on  the  satellite  S3-4.  The  UV  photometer  flown  on  the  same 
satellite  had  R  =  1200  counts  R'^  s'^.) 

Let  C  be  the  total  observed  count  rate,  D  the  observed  dark  count 
rate,  and  S  the  counting  rate  due  to  the  signal.  Then 

S  =  C  -  D  (1) 

If  the  instrument  counts  for  a  time  t,  then  the  statistical  uncertainties 


in  these  quantities  are 

Sc^  =  Ct 

(2) 

.  ot 

(3) 

$3^  =  Sj,^  +  Sp^  =  (C  +  D)t  =  (S  +  2D)t 

(4) 

Away  from  the  auroral  zone,  the  dark  count  rate  is  constant,  and  S  is 
simply  4itI  R,  where  R  is  the  previously  defined  instrument  response,  and 
47tl  is  the  apparent  column  emission  rate  of  the  signal. 

In  general,  we  wish  to  specify  a  maximum  acceptable  value 
relative  uncertainty,  u^  =  s^/S.  This  value  may  be  expressed 
R  and  t: 

for  the 

in  terms  of 
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(5) 


Solving  this  expression  for  R,  we  obtain  an  expression  for  the  instrument 
response  that  is  necessary  to  obtain  a  specified  statistical  uncertainty 


ir)  d  sp«i'ified  counting  time: 


1  +^1  +8u^^D1 
?(4r! )u^^t 


(6) 


This  expression  was  used  to  drive  the  values  given  in  Tables  4-6. 

Table  4  gives  required  R  values  for  various  integration  times, 
assuming  the  instrument  is  recording  the  minimum  1356  A  column  emission 
rates  listed  in  Table  3.  (0.1  Rayleighs  at  night  and  100  R  for  the 

daytime  and  aurora).  The  R  values  are  based  on  the  accumulation  of  25 
counts  of  real  signal  in  the  absence  of  either  dark  counts  or  background 
emission.  Thus,  Table  4  applies  to  ideal  conditions  leading  to  a 
statistical  uncertainty  of  20  percent.  An  uncertainty  greater  than  20 
percent  leads  to  unacceptable  uncertainty  in  the  deduced  electron 
dens i ty. 


Integration 
Time  (sec) 

Nighttime 
(0.1  Rayleighs) 

Day  and 

Aurora  (100  Rayleighs) 

.1 

2500 

2.5 

.2 

1250 

1.25 

.5 

500 

.5 

1. 

250 

.25 

2. 

125 

.125 

5 

50 

.05 

10 

25 

.025 

Table  4.  Instrument  Response,  R  (counts/R-sec ) ,  Required  to  Detect  the 
Indicated  Minimum  Expected  01  1356  A  Intensity  With  a  Statistical 
Uncertainty  of  20  Percent  or  Less.  (No  dark  counts.) 

Table  5  includes  required  R  values  similar  to  those  listed  in  Table  3 
but  for  more  realistic  conditions.  In  Table  5,  we  introduce  a  dark  count 
rate  of  6  count/s.  We  chose  this  value  based  on  examples  of  the  rate 
given  by  Huffman  et  al .  (1980).  Again,  the  given  R  values  lead  to  20 
percent  uncertainty. 
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Integration 
Time  (sec) 

Nighttime 
(0.1  Rayleighs) 

Day  and  Aurora 
(100  Rayleighs) 

.1 

2600 

2.61 

.2 

1360 

1.4 

.5 

600 

.60 

1 

340 

.34 

2 

200 

.20 

5 

no 

.11 

10 

70 

.069 

20 

45 

.045 

60 

25 

.025 

Table  5.  Similar  to  Table  2  Except  for  the  Presence  of  a  Dark  Count  Rate 
of  6  count/s. 


The  required  R  values  for  100  Rayleighs  of  01  1356  A  do  not  strain 
instrument  specifications,  even  for  short  Integration  times  like  0.1s. 

At  0.1  Rayleighs,  however,  the  necessary  counts  are  not  collected  until 
integration  times  approaching  one  minute  are  applied.  The  DMSP  satellite 
will  cover  a  distance  of  -  450  km  along  its  orbital  path  in  one  minute.  A 
single  data  point  over  450  Km  distance  will  be  unacceptable  in  situations 
where  is  noticeably  changing  on  shorter  distance  scales. 

The  R  values  in  Table  5  demonstrate  that  existing  UV  spectrometers 
cannot  reliably  detect  the  minimum  nighttime  01  1356  A  intensity. 
Therefore,  we  propose  the  inclusion  of  a  UV  photometer  similar  to  that 
flown  on  S3-4  (Huffman,  1980).  The  photometer  can  be  designed  either  to 
detect  only  01  1356  A  line,  or  to  detect  both  01  1356  A  and  01  1304  A 
simultaneously.  In  the  latter  case,  the  minimum  nighttime  intensity  of 
the  combined  features  is  about  0.3  R.  The  required  instrument  responses 
for  both  cases  and  for  various  integration  times  are  shown  in  Table  6. 

The  dark  count  rate  is  the  same  as  for  the  S3-4  photometer.  Since  the 
S3-4  photometer  had  R  =  1200  counts  R~^  s"^,  it  is  clear  that  existing 
technology  can  provide  accurate  measurements  with  reasonable  integration 
times.  Using  a  filter  wheel,  the  same  photometer  could  also  be  used  to 
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monitor  one  or  more  N2  LBH  bands.  These  bands  occasionally  obscure  01 
1356  A  at  night. 

Nighttime 


Integration  Dark  Count  Rate 


Time  (sec) 

No  Dark 

Counts 

of  4  counts/sec 

1356  alone 

1356  &  1304 

1356  alone 

1356  &  1304 

.05 

5000 

1700 

5100 

1700 

.1 

2500 

830 

2600 

860 

.2 

1300 

420 

1300 

440 

.5 

500 

170 

570 

190 

1 

250 

83 

310 

100 

2 

130 

42 

180 

60 

5 

50 

17 

93 

31 

10 

25 

8.3 

59 

20 

Table  6.  Instrument 

Response  Required  to  Detect  01  1356  A 

Alone  and  in 

Combination  with  01  1304  A. 


To  monitor  the  nighttime  6300  A  airglow,  we  propose  to  use  a  tilting 
filter  photometer  or  its  equivalent.  This  airglow  measurement  is 
complicated  by  the  presence  of  two  background  sources.  One  background 
source  is  OH  airglow  emission,  which  can  be  comparable  in  intensity  to  01 
6300  A.  The  other  source  is  reflected  light  from  astronomical  sources. 

At  new  moon,  this  source  is  almost  negligible.  At  full  moon,  reflected 
moonlight  can  be  as  much  as  100  times  larger  than  the  01  6300  A  line. 
Besides  subtracting  the  background,  we  must  also  correct  for  that  part  ot 
the  01  6300  A  intensity  that  is  reflected  from  the  surface.  Because  of 
these  complications,  we  must  discuss  the  instrument  design  in  more 
detai  1 . 

Because  the  astronomical  background  is  broadband,  it  is  most 
desirable  to  use  a  filter  with  the  most  narrow  possible  bandwidth. 

Filters  for  use  at  visible  wavelengths  can  be  made  with  bandwidths  of  2.0 
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to  ?.S  A.  Unfortunately,  thin  narrow  bandwidth  can  be  aciiieved  only  if 
the  field  of  view  is  about  0.01  sr  or  less.  To  subtract  the  background, 
without  introducing  an  unacceptable  statistical  uncertainty  and  with 
reasonable  integration  times,  we  may  have  to  increase  the  field  of  view 
by  about  an  order  of  magnitude.  Thus,  the  filter  bandwidth  will  be 
increased  to  5  to  10  A.  Introducting  of  optics  in  front  of  the  filter  to 
increase  the  effective  collecting  area  has  the  same  effect. 

Although  we  performed  an  extensive  quantitative  analysis  of  the 
detectability  of  01  6300  A  in  the  nighttime  ionosphere  under  this 
contract,  we  will  not  report  it  here  because  it  has  been  superseded  by  a 
more  recent  analysis  that  we  have  performed  under  a  contract  with  the 
Applied  Physics  Laboratory  (APL)  of  Johns  Jopkins  University.  The  latter 
analysis  is  based  on  more  complete  and  more  accurate  information 
concerning  the  01  63'''0  A  airglow,  the  sources  of  background,  and  optical 
insfrum--'nt  capabilities.  This  analysis  confirms  our  previous  conclusion 
tnat  thr-  measurement  is  feasible,  but  the  quantitative  results  are 
sufficiently  different  to  make  presenting  the  findings  of  the  old 
analysis  uninform  and  misleading.  The  new  analysis  will  appear  in  two 
reports  being  prepared  by  the  authors.  The  first  is  being  written  in 
collaboration  with  J.R.  Jasperse  of  AFGL  and  may  be  obtained  from  him 
when  it  is  complete.  The  second  is  the  final  report  for  the  APL  contract 
mentioned  above.  It  may  be  obtained  from  the  authors  when  it  is 
complete,  (In  both  cases  the  most  likely  availability  date  is  March 
1985. ) 

The  results  of  our  most  recent  analysis  may  be  summarized  as  follows: 

We  assumed  a  hypothetical  instrument  with  the  following  design 
parameters 

A  =  19.6  cm  (50  mm  filter,  no  optics) 

0.01  sr  (half  cone  angle  =  3°) 

=  bi; 

T  --  ?bk 

R  =  31?  cts  R'^  sec"^ 


l^w^  W  ■  ■  U  ■  t/»  ■  ■  IJ  '  <7"  w  "  r 


For  the  actual  local  time  of  the  "noon-midnight "  OMSP  orbit,  01  6300  A 
can  be  accurately  measured  by  our  hypothetical  instrument  under  all 
moonlight  conditions.  This  is  true  even  in  the  mid-latitude  trough.  The 
longest  integration  time  that  will  be  required  is  about  4.3  seconds  near 
40°  latitude  at  full  moon.  This  includes  the  time  required  for  observing 
the  background  intensities.  In  the  mid-latitude  trough,  the  longest 
integration  time  is  about  3.7  seconds.  Since  the  DMSP  satellite  moves  at 
about  7.5  km/sec,  this  implies  a  spatial  resolution  of  28  km.  As  the 
instrument  will  probably  have  a  field  of  about  0.01  sr,  the  actual 
spatial  resolution  will  be  about  85  km.  This  appears  to  be  quite 
adequate  to  map  the  density  gradients  of  the  trough  wall.  There  appears 
to  be  no  need  to  augment  the  photometer  with  a  telescope. 


:n 


5.  Preliminary  Analysis  of  Data  from  a  1978  Dayglow  Rocket  Experiment 

In  this  section,  we  describe  an  analysis  of  data  from  a  dayglow 
rocket  experiment  conducted  in  1978.  The  data  include  an  EDP  deduced 
from  ground-based  ionosonde  data  and  UV  intensities  measured  on  board  a 
sounding  rocket.  The  data  are  coincident  in  that  the  respective 
instruments  remotely  sensed  roughly  the  same  region.  This  analysis  came 
out  of  our  effort  to  find  experiments  conducted  in  the  dayglow, 
nightglow,  and  aurora  that  would  provide  both  EDP's  and  UV  intensities. 
The  purpose  for  the  search  was  to  test  our  predictive  capabilities,  which 
is  essential  to  the  concept  of  deducing  the  EDP  from  satellite  optical 
data.  We  must  have  confidence  that  our  models  can  be  used  to  predict  the 
EDP,  given  that  they  are  producing  good  optical  results.  We  have  found  a 
1978  dayglow  rocket  experiment  that  meets  our  criteria.  Comparison 
between  the  data  from  this  experiment  and  the  results  of  our  calculations 
has  revealed  that  we  cannot  reproduce  both  the  EDP  and  the  optical  data. 
We  do  not  have  the  answer  to  this  problem  which  is  most  pronounced  at  low 
altitudes  (<120  km).  The  same  problem  has  been  observed  in  similar 
comparisons  with  other  ionosonde  data.  We  plan  to  further  investigate 
the  problem  jointly  with  J.R.  Jasperse  and  D.W.  Decker  at  AFGL. 

The  rocket  experiment  has  been  discussed  by  Gentieu  et  al.  (1979), 
Meier  et  al .  (1980),  and  Anderson  et  al.  (1980).  An  analysis  of  the  data 
appears  in  the  Meier  paper.  The  rocket  was  launched  from  White  Sands 
Missile  Range  at  1300  MST  on  January  9,  1978.  The  solar  zenith  angle  was 
56°.  The  prime  experiments  involved  two  spectrometers  spanning  the 
region  from  530  to  1240  A.  Data  were  recorded  between  90  and  262  km  for 
nadir,  near  zenith,  and  tangential  viewing  directions.  The  ionosonde 
data  presented  in  this  section  were  given  to  us  by  Paul  Feldman  (private 
communication,  1983)  and,  to  our  knowledge  have  not  been  published. 

We  have  limited  our  calculations  to  the  EDP  and  associated  ion 
densities.  Our  optical  intensities  are  the  same  as  those  generated  by 
Meier  et  al .  (1980)  since  we  use  the  same  photoelectron  flux  code  and 
cross  sections  to  produce  the  altitude  profiles  of  volume  emission. 
Intensities  and  their  comparisons  with  data  for  selected  emissions  may  be 
seen  in  Meier  et  al.  (1980).  A  recent  analysis  by  Meier  and  colleagues 
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(unpublished)  has  refined  their  earlier  conclusions  regarding  atmospheric 
and  solar  conditions  that  were  present  during  the  1978  rocket  experiment. 
The  following  parameters  reflect  these  conditions  and  have  led  to 
excellent  agreement  between  predicted  intensities  and  the  spectrometer 
data  (01,  NI,  and  features): 

1)  960  K  Jacchia  (1977)  model  atmosphere 

2)  1.1  times  Torr  et  al .  (1979)  EUV  solar  spectrum  indexed  by  ^ 
=  71. 

3)  Donnelly  and  Pope  (1973)  soft  x-ray  solar  spectrum 

4)  0.6  times  Stone  and  Zipf  (1974)  1356  4  electron  impact  cross 
section 

We  have  used  the  information  in  parameters  1-3  to  generate  the 
calculated  EDP  shown  in  Figure  4.  The  calculations  were  done  with  our 
photoelectron  code  and  with  the  auroral  chemistry/optical  emission  code, 
which  was  modified  for  application  to  the  daytime  ionosphere.  A 
calculated  profile  by  Decker  and  Jasperse  (personal  communication)  is 
included  for  comparison.  They  applied  approximately  the  same  conditions 
as  listed  above.  Neither  profile  satisfactorily  fits  the  ionosonde  data, 
especially  at  the  lower  altitudes.  The  respective  results  came  from 
first  attempts  to  fit  the  data  and,  no  doubt  can  be  improved  with  certain 
parameter  adjustments.  Before  adjusting  parameters,  however,  we  have 
decided  to  resolve  the  discrepancies  between  the  two  models  producing 
these  results.  Resolving  the  discrepancies  is  being  done  by  comparing 
detailed  results  such  as  the  attenuated  solar  spectrum,  primary 
ionization  rates,  and  the  photoelectron  source  spectrum  at  selected 
altitudes  for  the  same  model  atmosphere  and  solar  spectrum.  A  simplified 
spectrum  is  being  used  to  facilitate  the  comparisons,  especially  those 
which  include  hand  calculations.  The  comparison  work  will  be  reported  at 
a  later  date.  We  can  say,  based  on  comparisons  to  date,  that  differences 
in  the  treatment  of  photoabsorption  and  photoionization  cross  sections 
are  likely  to  be  significant  in  producing  the  differences  shown  in  Figure 
4. 
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Figure  4.  Calculated  and  measured  EDP's  above  White  Sands  Missile  Range 
on  January  9,  1978  at  1300  MST.  The  calculated  results  are 
first  attempts  at  fitting  the  data. 
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In  suimary,  we  observe  that  our  predicted  EDP  lies  below  the  EDP 
deduced  from  the  ionosonde  data  by  a  factor  of  two  and  greater,  depending 
on  altitude.  At  the  same  time  our  predicted  intensities  for  features 
such  as  01  1356  A  and  LBH  bands  agree  with  the  observations.  To  address 
this  problem,  we  are  making  detailed  comparisons  with  the  results  of 
Decker  and  Jasperse  (personal  communication) and  will  be  examining  in 
detail  the  key  parameters  leading  from  ionization  rates  to  the 
corresponding  EDP. 
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Dependence  of  Auroral  Middle  UV  Emissions  on  the 
Incident  Electron  Spectrum  and  Neutral  Atmosphere 


Abstract 


In  this  paper  we  examine  the  relationship  between  two  important 
middle  UV  auroral  emissions,  the  electron  energy  spectrum,  and  the 
neutral  atmosphere.  Using  both  Maxwellian  and  Gaussian  forms  for  the 
incident  energy  spectra,  we  show  that  the  Vegard-Kaplan  band  system  is 
sensitive  to  both  the  atomic  oxygen  density  and  the  characteristic  energy 
of  the  incident  spectrum.  This  is  due  to  quenching  of  the  ^  state 

by  atomic  oxygen.  (Since  there  remains  some  uncertainty  about  the 
magnitude  of  the  quenching,  we  have  used  both  large  and  small  values  in 
our  calculations.)  We  also  show  that  the  second  po  itive  system  of  N2  is 
much  less  sensitive  to  incident  energy  and  atomic  oxyqen  density.  In 
principle,  observations  of  these  two  emissions  could  be  used  to 
characterize  the  incident  electron  spectrum  if  the  atomic  oxygen  density 
were  accurately  known  from  models  or  observations. 
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1.  Introduction 


In  this  paper,  we  present  results  showing  emission  properties  for 
selected  features  within  the  N2  Vegard-Kaplan  (VK)  and  second  positive 
(2P)  systems.  We  wish  to  know  how  such  features  vary  in  relation  to  one 
another,  and  to  other  features  as  well,  from  one  aurora  to  another.  The 
variation  results  from  changes  in  the  incident  electron  spectrum  and 
composition.  Our  interest  here  is  in  the  middle  ultraviolet  (MUV)  and 
for  this  region,  the  above  mentioned  systems  provide  the  most  prominent 
features,  with  the  exception  of  2972  A,  which  comes  from  0{^S).  This 
feature  will  not  be  discussed  here. 

Relative  variations  among  emission  features  offer  the  potential  to 
deduce  characteristics  of  the  energy  source  and  neutral  composition  from 
observations  of  their  intensities.  This  provided  the  motivation  for  the 
work  reported  by  Strickland  et  al .  [1983]  (hereafter,  SJW)  on  prominent 
far  UV  (FUV)  features.  The  work  being  reported  here  represents  an 
extension  of  that  work  into  the  MUV  region.  There  has  been  interest  over 
the  years  in  auroral  emission  variability,  mainly  as  observed  from  the 
ground.  Vallance  Jones  (e.g.  Vallance  Jones,  1975)  and  Rees  (e.g.  Rees 
and  Luckey,  1974)  have  pioneered  much  of  the  work  on  this  subject  (a 
brief  review  is  given  by  SJW).  Here,  we  have  in  mind  observations  from 
satellites  and  the  determination  of  their  usefulness  for  mapping  electron 
precipitation  characteristics  which  in  turn  provides  information  on 
electron  density  profiles  with  the  aid  of  models. 


2.  Atmospheric,  Molecular,  and  Chemical  Parameters 


We  have  used  the  same  model  atmosphere  as  was  used  by  SJW:  the  model 
of  Jacchia  [1977].  SJW  also  discussed  the  variability  of  neutral 
densities  and  compared  the  Jacchia  model  with  observations.  Following 
SJW,  we  have  examined  the  effects  of  reducing  the  0  density  by  as  much  as 
a  factor  of  three. 

The  features  we  have  investigated  are  the  Vegard-Kaplan  (VK)  and  the 
second  positive  (2P)  band  systems  of  N2.  A  simplified  energy  level 
diagram  for  the  triplet  states  of  is  given  in  Figure  1.  In 
calculating  Vegard-Kaplan  intensities,  we  have  made  extensive  use  of  the 
analysis  of  Cartwright  [1978].  He  found  that  the  low-lying  vibrational 
levels  of  the  A  state  were  populated  largely  by  cascading  from  the  B,  B', 
C,  and  W  states.  Because  cascading  takes  place  on  time  scales  that  are 
short  compared  to  the  lifetime  of  the  A  state  (about  2  seconds),  the 
production  of  any  of  these  states  is  tantamount  to  production  of  the  A 
state.  The  only  exception  is  the  C  state  which  predissociates 
approximately  half  of  the  time  [Porter  et  al.,  1976].  Therefore,  we  have 
calculated  the  total  production  of  all  levels  of  the  A  state  by  using  an 
effective  electron  impact  cross  section  a-j: 

"7  ‘  "a  ^°8  *  *  "W 

which  is  shown  in  Figure  2.  The  individual  cross  sections  were  taken 
from  Cartwright  et  al.  [1977].  We  apportioned  the  production  among  the  A 
state  vibrational  levels  according  to  the  population  rates  given  by 
Cartwright  [1978]. 
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In  particular,  we  assigned  2?  percent  of  the  production  to  the  v'=0  level 
and  16  percent  to  the  v'=l  level.  For  transition  probabilities,  we  have 
adopted  the  values  given  by  Sheiransky  [1969]. 

The  A  state  is  affected  by  chemistry  only  through  quenching.  We  have 

included  quenching  by  N2,  0^,  N,  0,  and  NO.  Unfortunately,  the  quenching 

rate  coefficient  for  the  most  important  species,  0,  is  also  the  most 

uncertain.  (See  Torr  and  Torr  [1982]  for  review.)  Sharp  [1971]  and  Torr 

-11  T  -1 

and  Sharp  [1979]  deduced  quenching  rates  of  9  x  10  cm  s  for  the 

-in  T 

v'=0  level  and  2  x  10  cm  s  ‘  for  the  v‘=l  level  (hereafter  called  the 
Sharp  values).  These  determinations  were  based  on  rocket  data. 
Ground-based  auroral  observations  by  Rees  et  al .  [1976]  and  Val lance 
Jones  and  Gattinger  [1976]  yielded  values  about  one-third  as  large. 
Laboratory  measurements  by  Piper  et  al .  [1981]  agreed  with  the  lower 
values:  2.8  x  10"^^  cm^  s"^  for  v'=0  and  3.4  x  10"^^  cm"^  s'^  for  v'=l 
(hereafter  called  the  Piper  values).  Recently,  McDade  and  Llewellyn 
[1984]  have  reanalyzed  several  sets  of  auroral  emission  data  using  the 
Piper  values.  Since  they  did  not  have  measurements  of  the  atomic  oxygen 
density,  they  had  to  use  a  model.  They  obtained  good  agreement  between 
theory  and  observation  using  the  CIRA  72  model  atmosphere,  which  provides 
support  for  the  Piper  values.  However,  the  fact  that  their  0  density 
came  from  a  model,  rather  than  from  observations,  prevents  the  question 
from  being  conclusively  settled.  Consequently,  the  calculations 
presented  in  the  next  section  were  performed  using  both  sets  of  values. 

The  2P  system  (arising  from  C  to  B  transitions)  is  easier  to  deal 
with.  Cascading  from  the  D  and  E  states  is  negligible,  so  the  production 
rate  consists  only  of  direct  excitation  by  electron  Impact.  The 
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transition  to  the  B  state  is  allowed,  so  quenching  is  unimportant.  The 
fraction  of  the  total  volume  emission  rate  which  appears  in  a  given  band 
is  obtained  from  the  product  q(v',0)  A(v',v'')  t(v')  where  q(v',0)  is  the 
Frank-Condon  factor  for  production  of  the  v'  level  of  the  C  state  from 
the  0  vibrational  level  of  the  ground  state,  A(v',v")  is  the  transition 
probability  for  the  v'-v"  transition,  and  t(v')  is  the  lifetime  of  the 
v'  level.  (We  assume  that  in  the  ground  state,  only  the  lowest 
vibrational  level  is  appreciably  populated.)  For  the  (0,0)  band,  this 
product  has  the  value  0.27  [Vallance  Jones,  1974].  The  predissociation 
mentioned  above  affects  only  the  higher  vibrational  levels  of  the  C 
state,  and  is  not  important  for  the  (0,0)  band  for  which  the  calculations 
of  the  next  section  were  done.  We  have  used  the  previously-cited  cross 
section  [Cartwright,  1977]  to  calculate  the  C  state  production  rate. 


38 


3.  Calculated  Volume  Production  Rates  and  Intensities 


As  in  SJW,  we  have  performed  calculations  for  both  Maxwellian  and 
Gaussian  incident  electron  energy  spectra  and  for  a  range  of 
characteristic  energies  typical  of  auroral  conditions. 


The  Maxwellian  spectrum  may  be  represented  by 


M 


where  F|^(E)  is  the  electron  flux  at  energy  E,  is  the  total  energy 
flux,  and  is  the  characteristic  energy.  We  assume  that  the  flux  is 
isotropic  over  the  downward  hemisphere. 

The  Gaussian  spectrum  may  be  represented  by 


where  Fg,  and  Qg,  and  Eg  are  analogous  to  and  E^.  W  is  the  width 

of  the  spectrum.  (To  be  exact,  the  full  width  at  half  maximum  is 
2Vln  2  W.) 


Note  that  the  mean  energy  of  a  Maxwellian  spectrum  is  twice  the 
characteristic  energy:  <E>j,,  =  2  Ej,|.  In  contrast,  the  mean  energy  of  a 
Gaussian  spectrum  is  equal  to  the  characteristic  energy.  <E>g  *  Eg.  This 
should  be  kept  In  mind  when  comparing  calculations  based  on  the  two 


different  spectral  forms.  In  all  calculations  presented  here,  we  have 

?  1 

chosen  and  to  be  1  erg  cm  s  .  Since  volume  emission  rates  are 
directly  proportional  to  Q,  it  is  easy  to  scale  the  results  to  other 
values  of  the  energy  flux. 

The  altitude  dependence  of  the  production  and  loss  terms  for  the  lowest 
vibrational  level  of  N2(A)  are  shown  in  Figures  3  and  4.  The  Piper 
values  for  quenching  by  atomic  oxygen  were  used  in  Figure  3  and  the  Sharp 
values  were  used  in  Figure  4.  The  main  effect  of  the  higher  Sharp  values 
is  to  reduce  the  peak  emission  rate  by  a  factor  of  two.  A  less 
significant  effect  is  the  lifting  of  the  altitude  of  peak  emission  by 
about  5  kilometers.  Altitude  profiles  of  the  volume  emission  rate  of  the 
(0,6)  Vegard-Kaplan  band  at  2762  A  for  various  values  of  the 
characteristic  electron  energy  are  given  in  Figures  5  and  6.  Note  that 
even  though  the  peak  production  rate  increases  with  increasing  electron 
energy  (see  Figure  3  and  4),  quenching  causes  the  peak  volume  emission 
rate  to  decrease. 

Because  quenching  is  not  important  for  the  C  state,  the  volume 
emission  rate  for  any  band  of  the  2nd  positive  system  is  proportional  to 
the  production  rate.  In  particular,  the  volume  emission  rate  of  the 
(0,0)  band,  which  is  not  affected  by  predissociation,  is  equal  to  the 
production  rate.  The  altitude  dependence  of  the  (0,0)  band  volume 
emission  is  shown  is  Figure  7.  Note  that  in  contrast  to  the 
Vegard-Kaplan  system,  higher  energy  fluxes  produce  greater  peak  volume 
emission. 
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A.  Calculated  Column  Emission  Rates 


We  have  integrated  the  volume  emission  rates  shown  in  Figures  5-7  to 
proauce  column  emission  rates  as  functions  of  electron  energy.  The 
results  for  the  Vegard-Kapl an  system  are  shown  in  Figures  8  and  9.  We 
show  two  of  the  brightest  bands  of  the  system,  as  well  as  the  sums  of  the 
v'=0  and  v'=l  progressions.  Except  at  the  lowest  energies,  the  intensity 
of  the  Vegard-Kaplan  system  decreases  rapidly  with  increasing  electron 
energy.  This  is  the  effect  of  quenching.  Since  the  quenching  rate  is 
proportional  to  the  atomic  oxygen  concentration,  a  reduction  in  the 
quenching  coefficient  by  any  factor  is  equivalent  to  an  increase  in 
atomic  oxygen  by  the  same  factor.  Thus,  these  two  figures  can  also  be 
taken  as  illustrations  of  the  sensitivity  of  the  emission  to  variations 
in  oxygen  density.  That  is.  Figure  9  can  be  taken  as  the  intensities 
that  would  be  obtained  if  the  Piper  values  were  used,  but  with  the  oxygen 
density  increased  by  a  factor  of  three.  (This  neglects  the  second  order 
effects  due  to  the  additional  absorption  of  energy  by  the  atomic  oxygen, 
as  it  is  increased  relative  to  that  of  N^. ) 

The  energy  dependence  of  the  column  emission  rates  for  the  (0,0)  band 
of  the  2nd  positive  system  is  shown  in  Figure  10.  Although  this 
intensity  is  a  monotonically  increasing  function  of  energy,  the  rate  of 
increase  is  not  very  rapid,  especially  at  higher  energies.  At 
sufficiently  low  electron  energies,  most  of  the  energy  is  deposited  at 
high  altitudes  where  atomic  oxygen  dominates.  As  the  energy  increases,  a 
greater  fraction  is  deposited  in  nitrogen,  so  the  column  emission  also 
increases.  However,  when  the  electron  energy  is  high  enough,  the  energy 
is  deposited  at  altitudes  where  atomic  oxygen  is  not  important. 


so  the  column  emission  depends  mainly  on  the  total  energy  flux,  and  not 

on  the  characteristic  energy.  Although  the  peak  volume  emission 

rate  does  Increase  with  characteristic  energy,  the  emitting  layer  becomes 


thinner. 


The  two  band  systems  discussed  in  this  paper  could  not  be  used  to 
deduce  the  incident  electron  spectrum  from  optical  observations  without 
additional  information.  If  the  atomic  oxygen  density  were  known  from 
other  observations  or  from  models,  then  these  two  emissions  could,  in 
principle,  be  used  to  calculate  the  energy  content  and  mean  (or 
characteristic)  energy  of  the  electron  spectrum.  This  is  because  the  two 
emissions  have  different  dependences  on  characteristic  electron  energy, 
and  are  both  proportional  to  the  total  energy  flux.  The  main  problem  in 
making  use  of  these  emissions  is  their  relatively  low  intensities. 
However,  there  are  some  other  complications.  One  of  these  is  the  finite 
width  of  the  bands.  The  rotational  lines  making  up  a  band  are  spread 
over  several  Angstroms.  The  amount  of  spread  depends  on  the  rotational 
temperature  of  the  emitting  species.  For  example,  if  the  effective 
rotational  temperature  of  in  the  A  state  is  800  K,  the  VK  (0,6)  band 
is  spread  over  approximately  60  A.  (  We  have  not  attempted  to  model 
rotational  distributions  for  actual  auroral  conditions).  Rotational 
spreading  also  causes  the  VK  (0,9)  band  to  overlap  the  2P  (0,0)  band. 
Conway  and  Christensen  [1984]  modeled  this  contamination  of  photometer 
observations  in  the  dayglow  and  concluded  that  2P  (0,0)  band  observations 
should  be  done  spectrometrically  rather  than  photometrically.  This 
conclusion  is  almost  certainly  valid  for  auroral  conditions  as  well. 

There  is  a  second  complication  for  the  (0,0)  band  of  the  2nd  positive 
system:  Rayleigh  scattering  from  the  atmosphere  below  the  emitting 
layer.  This  is  not  a  problem  at  wavelengths  shorter  than  3000  A  because 
of  pure  absorption  by  Op-  wavtlength  of  the  (0.0)  band  (3371  A). 


43 


there  is  very  little  absorption,  but  there  is  sufficient  multiple 
scattering  that  most  of  the  light  emitted  in  the  downward  direction  is 
scattered  upward.  Thus  the  apparent  column  emission  rate  inferred  from 
satellite  observations  will  be  too  large  unless  this  effect  is  taken  into 
account.  (These  statements  are  based  on  multiple  scattering  calculations 
of  Anderson  quoted  in  Strickland  and  Oaniell  [1982].) 
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Figure  Captions 


Figure  1.  A  simplified  energy  level  diagram  for  the  triplet  states 
and  band  systems  of  molecular  nitrogen.  The  two  band  systems 
discussed  in  this  paper  are  the  Vegard-Kaplan  and  2nd  positive 
sy terns . 


Figure  2.  The  total  cross  section  for  the  production  of  the  A 

state  of  Nz  by  electron  impact.  This  cross  section  takes  into 
account  cascading  from  higher  lying  states.  (See  text  for 
complete  explanation.) 


Figure  3.  Altitude  profiles  of  the  production  and  loss  rates  for  the 
v'  =  0  vibrational  state  of  the  A  state  of  Na  for  (a)  Maxwellian 
and  (b)  Gaussian  incident  electron  spectra.  The  quenching  rate 
was  calculated  using  the  Piper  et  al .  [1961]  values  for  the 
quenching  coefficients. 


Figure  ^t.  Altitude  profiles  of  the  production  and  loss  rates  for  the 
v'  -  0  vibrational  state  of  the  A  state  of  Na  for  (a)  Maxwellian 
and  (b)  Gaussian  incident  electron  spectra.  The  quenching  rate 
was  calculated  using  the  Sharp  [1971]  and  Torr  a:.J  Sharp  [1979] 
values  for  the  quenching  coefficients. 


Figure  5.  Altitude  profiles  of  the  volume  emission  rate  of  the  (0,6) 
Vegard-Kaplan  system  for  (a)  Maxwellian  and  (b)  Gaussian  incident 
electron  spectra.  Profiles  are  shown  for  several  characteristic 
energies.  The  Piper  quenching  coefficients  were  used. 
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Figure  6.  Altitude  profiles  of  the  volume  emission  rate  of  the  (0,6) 
Vegard-Kaplan  system  for  (a)  Maxwellian  and  (b)  Gaussian  inct  ent 
electron  spectra.  Profiles  are  shown  for  several  characteristic 
energies.  The  Sharp  quenching  coefficients  were  used. 


Figure  7.  Altitude  profiles  of  the  volume  emission  rate  of  the  (0,0) 
band  of  the  2nd  positive  system  for  (a)  Maxwellian  and  (b) 
Gaussian  Incident  electron  spectra.  Profiles  for  several 
characteristic  electron  energies  are  shown. 


Figure  8.  The  dependence  of  the  column  emission  rate  of  the  (0,6) 
Vegard-Kaplan  band  on  the  incident  electron  energy  for 
(a)  Maxwellian  and  (b)  Gaussian  electron  spectra.  The  Piper 
quenching  coefficients  were  used. 


Figure  9.  The  dependence  of  the  column  emission  rate  of  the  (0,6) 
Vegard-Kaplan  band  on  the  incident  electron  energy  for 
(a)  Maxwellian  and  (b)  Gaussian  electron  spectra.  The  Sharp 
quenching  coefficients  were  used. 


Figure  10.  The  dependence  of  the  column  emission  rate  of  the  (0,0) 
2nd  positive  band  on  the  incident  electron  energy  for  (a) 
Gaussian  and  (b)  Maxwellian  electron  spectra. 
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This  paper  addresses  (he  probtea  at  monitoring  daytime  and  auroral  ionospheres  using  optical  measurements 
from  satellites.  The  most  attractire  approach  is  to  use  optical  emission  leatures  o(  the  neutral  atmosphere  to 
monitor  both  some  variability  and  neutral  density  and  composition.  From  this  information,  the  electron  density 
profile  may  be  deduced  subject  to  uncertainties  caused  by  neutral  winds  and  electric  HeMs.  Emission  features, 
which  we  examine  as  candidates  for  use  in  iMs  approach,  are  bands  of  the  N2  LBH  system,  N2  Vefjjrd- Kaplan 
2762  K,  fif  IN  3414  A,  and  the  atomic  oxygen  lines  at  I3S6  and  2972  A.  (.'alculated  intensities  of  these  features 
and  <1,  arc  shown  for  daytime  and  auroral  conditiotu. 


Introduction 

N  ttiis  paper,  we  shall  discuss  work  in  progress,  the  goal  of 
which  is  to  assess  the  utility  of  satellite-observed  optical 
emission  features  in  deducing  the  E  and  F  region  electron 
density  distributions  (n,).  Here  we  restrict  ourselves  to  the 
daytime  E  and  F  regions  and  to  the  auroral  E  region.  We 
know  of  no  emission  feature  in  these  regions  that  ^ves  a 
direct  signature  of  the  ions  present,  unlike,  say,  the  O  ‘  recom¬ 
bination  emissions  in  the  tropical  nighttime  airglow.'-^  O*  834 
A  comes  closest  to  fulfilling  the  requLrement  since  can 
affect  the  834  A  intensity  through  multiple  scattering.’  This  in 
itself  poses  a  difficult  problem  which  we  shall  not  address 
here.  It  then  becomes  a  two-step  process  to  obtain  n,  from 
optical  measurements.  We  must  first  obtain  a  representation 
of  the  source  (solar  ionizing  spectrum  or  incident  auroral 
electron  spectrum),  followed  by  its  use  in  calculating  n,.  The 
emphasis  in  this  paper  will  be  on  how  well  the  source  spec¬ 
trum  can  be  determined  from  optical  data.  To  address  the 
second  step,  we  give  examples  of  the  sensitivity  of  n,  to 
variations  in  the  source  spectrum. 

The  variability  of  n,  can  pose  some  serious  problems  for 
the  above  technique.  This  variabtUty  is  caused  by  1)  time  and 
spatial  changes  in  the  source,  2)  the  same  type  of  changes  in 
the  major  neutral  density  profiles,  and  3)  mechanics  which 
produce  bulk  transport  of  the  plasma.  As  noted  above,  item  1 
shall  receive  most  of  our  attention.  Item  2  is  important  be¬ 
cause  uncertainties  introduced  can  affect  both  the  accuracy  of 
the  deduced  source  spectrum  and  the  subsequently  calculated 
n,  profile.  The  accuracy  of  n,  is  dependent  not  only  on  the 
assumed  source  spectrum,  but  also  on  assumed  neutral  densi¬ 
ties  through  chemistry  and  diffiLsion.  More  shall  be  said  about 
item  2  in  the  third  section.  Item  3  has  to  do  with  winds  and 
fields  and  will  be  briefly  discussed  in  the  next  section.  Beyond 
this,  however,  we  shall  restrict  ourselves  to  conditions  for 
which  n,  is  dependent  only  on  the  ionizing  source. 

We  shall  examine  N2  LBH  1383  A  and  OI  1356  A  in  the 
dayglow.  For  auroral  conditions,  we  shall  consider  these  as 
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well  as  01  2972  A;  the  LBH  btmds  at  1273, 1325,  and  1354  A; 
Nj  Vegard-Kaplan  (VK)  2762  A;  and  N2*  IN  3914  A.  Some  of 
the  auroral  results  to  follow  have  recently  been  discussed  in 
more  detail  in  the  paper  by  Strickland  et  al.^  In  particular, 
these  authors  determined  the  variability  of  nadir  viewing 
sateUite-observed  intensities  for  OI  1356_,  3914  A,  and 

LBH  bands  at  1273,  1325,  and  1383  A  as  a  function  of 
electron  spectral  hardness  and  model  atmosphere.  The  causes 
of  variability  in  the  re'  live  intensities  among  the  above  fea¬ 
tures  are  variations  in  1 )  excitation  efficiencies,  2)  pure  photon 
absorption  (here  due  to  Oj),  and  3)  chemied  processes 
(quendiing,  in  particular)  that  are  due  to  changes  in  source 
and  neutral  atmospheric  conditions.  Results  will  follow  that 
demonstrate  the  strengths  of  these  effects  on  the  emissions. 

Ionospheric  Variability  Due  to  Neutral 
Winds  and  Electric  Fields 

In  the  F  region  ionosphere  (above  180  km)  both  the  ion- 
neutral  and  electron-neutral  collision  frequencies  are  much 
less  than  the  respective  gyrofrequcncies.  This  means  that  the 
mobility  of  a  charged  particle  is  much  greater  along  a  field  line 
than  across  field  lines.  The  neutral  winds  of  the  thermospheric 
circulation  produce  a  frictional  force  on  the  ionospheric 
plasma,  but  only  the  component  of  that  force  which  acts  along 
a  field  line  is  effective  at  imparting  motion  to  the  plasma. 
Magnetic  field  lines  lie  approximately  in  meridional  planes 
and  rise  in  altitude  toward  the  equator.  During  daytime  hours, 
the  normal  thermospheric  circulation  is  from  equator  to  pole,’'* 
which  tends  to  drive  the  plasma  downward.  This  reduces  not 
only  the  altitude  of  the  ionization  peak,  but  also  the  peak 
concentration,  since  chemical  losses  are  larger  at  lower  alti¬ 
tudes.  At  middle  latitudes  the  effects  of  these  winds  can  be  as 
important  as  diffusion,  resulting  in  a  lowering  of  the  F2  peak 
by  -  50  km  and  a  reduction  in  peak  electron  density  by 
-  35^1.’'*  However,  the  thermospheric  circulation  is  quite 
variable  on  a  daily  basis,  which  makes  the  ionosphenc  varia¬ 
bility  difficult  to  predict.  One  source  of  variability  is  the 
energy  deposited  at  high  latitudes  during  geomagnetic  sub¬ 
storms  which  can  alter  the  global  thermospheric  circulation 
and  even  reverse  the  direction  of  flow.*-’  Furthermore,  the 
resistance  of  ions  to  cross-field  motion  can  alter  the  thermo¬ 
spheric  winds  in  the  altitude  regime  above  300  km,  further 
compUcating  the  picture. 

In  the  E  region  iono.sphere  (150-160  km)  the  ion  cross-field 
mobility  is  considerably  larger  than  that  of  electrons.  In  this 
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altitude  regime,  neutral  winds  induce  a  charge  separation 
which  results  in  a  global  electric  field  system.  Because  mag¬ 
netic  held  lines  are  nearly  equipoienhals.  this  electric  held 
system  is  map^d  into  the  F  region,  where  it  causes  plasma 
{Ex  B)  drifts.^  One  may  think  of  the  E  region  winds  acting 
as  a  dynamo  generator  while  the  F  region  plasma  responds  as 
an  electric  motor.  The  effect  of  the  electric  held  is  generally 
smaller  than  that  due  to  F  region  winds  and  varies  throughout 
the  day. 

In  addition  to  the  large-scale  effects  caused  by  thermo¬ 
spheric  winds  and  the  dynamo  electric  held,  there  are  small- 
scale  structures— ionospheric  irregularities— which  have  a 
variety  of  sources.’  The  irregularities  are  present  in  both  the  E 
and  F  regions  and  at  all  latitudes."’"  At  auroral  and  equa¬ 
torial  latitudes  the  primary  causes  of  ionospheric  irregularities 
are  plasma  instabilities,  although  the  relative  importance  of 
neutral  dynamics  in  the  explanation  of  equatorial  spread  F  is 
subject  to  debate.  At  midlatitudes,  the  situation  is  less  well 
understood.  Wind  shears  are  known  to  produce  sporadic  E 
layers  while  traveling  ionospheric  disturbances  have  been  ex¬ 
plained  by  gravity  waves  propagating  from  below  or  from  the 
auroral  zone.  A  number  of  plasma  mechanisms  have  been 
proposed  to  explain  smaller-scale  structures,  but  veriheation 
of  these  hypotheses  awaits  more  complete  sets  of  ob.serva- 
tional  data  " 

Both  E  and  F  region  thermospheric  winds  can  have  a 
significant  influence  on  the  c  Icciron  distribution  above  200  km 
or  so.  They  are  hard  to  observe  and  continue  to  be  the  subject 
of  theoretical  and  experimental  research  If  the  winds  are 
known,  their  effects  on  the  ionosphere  can  be  determined/ 
The  principal  implication  for  satellite  monitoring  of  the  iono¬ 
sphere  is  that  some  method  of  measuring  or  calculating  ther¬ 
mospheric  winds  must  be  developed  Global  models  of  the 
thermospheric  circulation  provide  a  useful  beginning,  but  the 
effects  of  geomagnetic  activity  must  al.so  be  included  if  rea¬ 
sonable  accuracy  is  to  be  obtained.  Since  ionospheric  irregu¬ 
larities  are  less  well  understood,  particularly  at  midlatitudes, 
they  present  a  formidahle  problem  for  satellite  monitoring. 
Those  problems  are  not  addressed  further  in  this  paper,  but 
remain  important  research  subjects. 

(Emission  Variability  due  to  Sources 
and  Neutral  Atmosphere 

The  technique  being  addressed  is  most  attractive  for  situa¬ 
tions  in  which  the  den.sity  profiles  of  N^,  O2,  and  O  are 
known  Under  those  conditions,  variability  of  one  optical 
intensity  relative  to  another  with  changes  in  either  time  or 
location  can  be  directly  related  to  variability  in  source  condi¬ 
tions  However,  the  neutral  density  profiles  cannot  be  pre- 
ci.scly  known  because  of  their  intrinsic  variability  and  the 
difficulty  of  precisely  descnbing  that  vanability.  The  latter 
problem  is  due  to  the  approximate  nature  of  existing  thermo¬ 
spheric  den.sity  models  (for  recent  models,  sec  Refs.  12-16) 
and  to  the  limitations  of  measurement  techniques.  Of  particu¬ 
lar  interest  here  arc  satellite-based  measurements.  In  situ 
measurements  can  he  expected  to  be  the  most  reliable,  but 
may  be  insufficient  for  .specifying  altitude  profiles  to  the 
degree  of  accuracy  desired  Optical  remote  sensing  techniques, 
on  the  other  hand,  do  hold  a  potential  for  measuring  column 
den.sitics,  which  can  be  used  to  scale  corresponding  yolume 
density  profiles  This  problem  is  addre.s.sed  next.  One  of  the 
objectives  of  the  present  work  must  be  to  determine  just  how 
precisely  the  major  neutral  densities  need  to  be  known  to 
make  the  technique  iindei  investigation  practical  The  required 
precision  is  dependent  on  .source  conditions  and,  of  course,  on 
instrumental  conditions 

The  problem  raised  above  is  whether  it  is  possible  to 
separate  the  efleets  of  neutral  density  variations  from  source 
variations  when  using  optical  data  Various  re.searchers  are 
actively  investigating  the  problem  of  sensing  the  neutral  atmo¬ 
sphere  '  Their  roncern  has  been  with  the  analysis  of  satel¬ 


lite-observed  limb  profiles  for  features  such  as  01  1356  A  and 
Nj  LBH  bands.  We  have  concentrated  our  efforts  on  nadir 
viewing  observations,  and  we  bcheve  that,  with  the  right 
choice  of  features  and  instrumental  parameters,  the  source 
and  atmospheric  variabilities  can  Ik  decoupled.  This  is 
strengthened  by  the  fact  that  there  are  limitations  to  the 
relative  variability  between  the  N2  and  O2  densities  in  the 
regions  of  peak  excitation  (for  either  auroral  or  daytime 
conditions).  Furthermore,  the  models  dted  above  are  probably 
more  than  adequate  under  most  observing  conditions  for 
describing  these  relative  densities.  It  then  becomes  attractive 
to  use  selected  N2  band  emissions  for  initially  characterizing 
the  source  spectrum.  For  auroral  conditions,  their  relative 
strengths  indicate  the  hardness  of  the  electron  source  spec¬ 
trum  through  cither  of  the  basic  mechanisms  of  pure  absorp¬ 
tion  or  quenching,  depending  on  the  features.  Given  a  repre¬ 
sentation  of  the  spectral  hardness,  the  magnitudes  of  the 
observed  intensities  then  determine  the  energy  content  of  the 
spectrum.  In  the  daytime  ionosphere,  most  of  the  F  region 
ionization  is  produced  by  radiation  in  the  range  l(X)-800  A 
with  a  smaller  contribution  from  the  8(X)-910  A  range.*  The 
radiation  below  ISO  A  produces  only  5  to  8%  of  the  ioniza¬ 
tion,  depending  on  solar  activity.  Since  the  shape  of  the  solar 
spectrum  in  the  wavelength  range  150-910  A  does  not  vary 
drastically  with  solar  activity,"  the  important  parameter  is  the 
energy  content  of  that  part  of  the  spectrum.  The  magnitudes 
of  the  observed  airglow  intensities  provide  a  measure  of  that 
parameter.  However,  the  variability  of  the  ultraviolet  dayglow 
is  not  well  understood.  Rocket  data  analyzed  by  Meier  et  ai. 
(Ref.  20  and  personal  communication)  indicate  that  the  varia¬ 
bility  is  weaker  than  the  results  presented  below  would  sug¬ 
gest.  This  remains  an  area  of  active  investigation. 

With  an  initial  characterization  of  the  source  spectrum,  we 
may  turn  to  the  problem  of  better  determining  the  neutral 
densities.  The  atomic  oxygen  density  is  of  prime  interest,  since 
it  probably  experiences  the  greatest  relative  variation  (see 
above  papers  on  thermospheric  modeling)  and  has  associated 
wth  it  some  key  uv  emission  features  (examples  are  01  1 304 
A,  01  1356  A,  and  01  2972  A).  Because  O  is  a  minor  species 
when  considered  over  the  entire  excitation  region,  its  optical 
intensities  are  sensitive  to  its  overall  density.  Thus,  given  the 
source  spectrum,  a  measure  of,  say,  the  01  1356  A  intensity 
(with  a  minor  contamination  from  N2  LBH  1354  A)  gives  a 
direct  measure  of  the  magnitude  of  the  O  density  profile.  We 
may  now  repeat  the  process  to  further  improve  both  the 
source  representation  and  neutral  density  description. 

Computational  Models 

Daytime  Ionosphere 

Wc  possess  two  models  and  associated  computer  codes  for 
describing  the  interactions  of  solar  ionizing  radiation  and 
photoelectrons  with  the  aUnosphere.  One  is  the  Boltzmann- 
Fokker-Planck  model  of  Jasperse,^'-’’  which  provides  de¬ 
tailed  energy  spectra  of  the  electrons  over  both  the  thermal 
and  photoelcctron  energy  regions.  The  other  model  is  by 
Strickland  and  Meier, and  also  gives  photoelcctron  spectra, 
but  with  a  less  detailed  description  of  the  Boltzmann  collision 
integral.  In  the  second  model,  more  emphasis  is  placed  on  the 
calculation  of  optical  intensities  Chemistry  codes  arc  coupled 
to  the  abovementioned  photoelcctron  codes  to  calculate  tem¬ 
peratures,  ion  densities,  and  various  neutral  densitie.s  Results 
to  follow  will  come  from  the  use  of  the  second  photoelcctron 
model  mentioned  above  and  associated  chemistry  code  which 
allows  for  diffu.sion 

Auroral  Ionosphere 

The  electron  transport  model  of  Strickland  el  al  is  used 
to  speafy  the  electron  spectrum  and  a.s.sociatcd  excitation  and 
ionization  rates.  A  Boltzmann  equation  is  solved  which  gives 
the  spectrum  or  distribution  function  from  the  highest  primary 
energy  of  interest  down  to  the  eV  range.  EJectron  interactions 
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are  limited  to  the  particle- particle  type,  which  iocludes  elastic 
scattering,  excitation,  and  ionization  leading  to  energy  loss  of 
the  incident  electron  and  production  of  a  secondary  electron. 
A  time-dependent  chemistry  code  is  coupled  to  the  transport 
re.suJts  and  gives  the  densities  of  numerous  ions  and  neutrals, 
along  with  intensities  for  a  variety  of  emission  features.  This 
code  is  essentially  the  same  as  the  one  from  which  the  daytime 
results  to  follow  were  obtained. 

Comparisons  with  Optical,  Electron,  and  Ion  Data 

The  calculation  of  photoelectron  and  auroral  electron  spec¬ 
tra,  excitation  rates,  chemical  species,  densities,  and  optical 
intensities  involves  a  large  body  of  input  data  and  the  use  of 
some  large  computer  codes.  It  is  essential  to  compare  these 
calculated  quantities  with  independent  results  as  well  as  mea¬ 
surements  prior  to  a  serious  investigation  of  the  problem  being 
addressed.  Three  such  comparisons  follow. 

We  first  note  that  electron  spectra  are  the  most  difficult  to 
calculate  the  abovemeniioned  quantities.  We  have  tested  them 
through  various  comparisons  and  by  observing  how  well  they 
conserve  energy.  The  latter  test  typically  yields  an  error  of  less 
than  10%  The-auroral  electron  results  have  been  well  tested, 
beginning  with  comparisons  discussed  by  Strickland  et  al.^ 
followed  by  subsequent  comparisons  with  rocket  and  satellite 
data.  This  latter  work  is  unpublished  but  does  generally  show 
good  agreement.  Specific  examples  of  applied  data  are  the 
rocket  data  of  Spiger  and  Anderson'**  and  satelhte  data  of 
Craven  and  Frank.**  The  observed  quantity  is  the  auroral 
electron  spectrum  given  as  a  function  of  Irath  energy  and 
angle.  With  regard  to  the  calculated  photoelectron  spectra, 
excellent  agreement  has  been  achieved  with  both  data  and 
independent  calculations.  These  comparisons  may  be  seen  in 
the  report  of  Strickland  and  Meter.** 

For  the  comparisons  to  follow,  we  have  chosen  optical,  ion 
density,  and  electron  density  data,  since  these  are  the  types  of 
data  we  are  ultimately  interested  in.  We  begin  with  the  elec¬ 
tron  density  information  in  Fig.  1,  which  comes  from 
Jasperse.  **  This  figure  shows  good  overall  agreement  between 
theory  and  experiment  for  daytime  conditions.  An  example  of 
fractional  ion  abundances  under  auroral  conditions  is  shown 
in  Fig.  2.  The  data  come  from  Swider  and  Narcisi*'  and  were 
obtained  with  an  ion  mass  spectrometer  mounted  on  a  rocket 
fired  into  the  continuous  aurora  The  calculated  results  show 
good  overall  agreement,  with  the  exception  of  those  for  N2* 
which  lie  above  the  data  For  this  same  experimenL  satisfac- 
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Fig.  2  Fractional  ion  abundances  lor  the  continuous  aurora.  Data  are 
given  by  (he  solid  and  broken  curves  and  come  from  Swider  and 
Narcisi.  “ 
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Rg.  3  1356  A  dayglow  limb  intensity.  Data  are  from  Newman  and 
Clrtslensen”’'*  while  calculations  are  from  Strickland  and  Anderson.* 


tory  agreement  was  also  obtained  between  the  calculated 
electron  density  and  coincident  radar  data.  As  an  example  of 
a  comparison  with  optical  data,  we  show  in  Fig.  3  results  from 
Strickland  and  Anderson*’  which  give  the  1356  A  limb  inten¬ 
sity  for  daytime  conditions.  The  data  were  provided  by 
Newman  and  Giristensen  from  a  spectrometer  experiment  on 
board  satellite  DMSP-F4,‘*'’  Calculations  of  the  needed 
photoclectron  spectra  and  excitation  rates  come  from  the 
second  of  the  previously  discus.sed  photoeicctron  models,** 
while  the  remaining  calculatiotis  were  provided  by  the  photon 
transport  model  of  Anderson  et  al. 

Daytime  Results 

In  this  section,  we  wish  to  demonstrate  the  vanability  of 
optical  emissions  and  n,  caased  bv  changes  in  the  solar 
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Fig.  4  Zenith  and  nadir  viewing  01  1356  A  and  N]  LBH  1383  A 
dayglow  intensitie^.  The  solar  aenith  angle  is  60  deg.  The  labeling 
refers  to  the  following  conditions;  1)  low  solar  activity,  Jacchia  1977 
11(0);  2)  high  solar  activity,  Jacchia  1977  ii(0);  3)  low  solw  activity,  half 
of  Jacchia  1977  «(0). 
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fig.  5  Cdculated  di>1iinc  electron  density.  Labeling  is  tbe  same  as  in 
Fig.  4. 


activity  and  the  neutral  dcn.sities  We  have  ju.st  begun  addrcs.s- 
ing  this  problem  and  will  restrict  the  results  to  two  key  optical 
features  and  three  starting  conditions.  The  features  to  be 
considered  are  LBH  1.383  A  (the  (2,0)  and  (5,2)  bands)  and 
Ol  1 356  A  Both  features  are  at  sufficiently  short  wavelengths 
that  Rayleigh  scattering  of  solar  radiation  does  not  oon- 
tainmaic  their  intensities  on  emerging  from  the  ionosphere. 
Conditions  are  defined  in  terms  of  solar  activity  (indexed  by 
F,n,,  the  10  7  cm  .solar  flux)  and  the  model  aU7io.sphcre  We 
use  solar  spectrum  representations  from  Ton  et  al  for  high 
(F|„,-  207)  and  low  (Fj||,«7i)  activity.  The  basic  model 
atmosphere  comes  from  the  1977  Jacchia  model  One  set  of 
results  will  be  shown  for  which  the  O  density  has  been 
reduced  by  half 


Figure  4  shows  nadir  and  zenith  viewing  intensities  for  the 
three  sets  of  starting  conditions  We  observe  a  change  in 
intensity  of  ~  2  in  going  from  low  to  high  solar  activity.  As 
expected,  halving  the  atomic  oxygen  density  reduces  the  Ol 
1356  A  intensity  by  a  similar  factor.  The  LBH  1383  A  inten¬ 
sity  is  altered  only  slightly  because  of  a  small  reapportionment 
of  photoelcctron  energy  received  by  Nj.  As  we  mentioned 
^ve,  the  actual  variation  in  Ol  1356  A  due  to  solar  variabil¬ 
ity  may  be  less  than  our  calculations  indicate.  If  so,  it  may  be 
that  the  Torr  et  al.’*  solar  spectra  overestimate  the  amount  of 
euv  variabihty.  This  problem  remains  under  investigation. 

Figure  5  shows  electron  densities  for  the  two  levels  of  solar 
activity.  The  results  above  200  km  include  the  effect  of  diffu¬ 
sion  We  regard  these  results  as  preliminary,  since  they  are  the 
first  we  have  obtained  for  daytime  conditions  and  since  the 
diffusion  model  is  still  under  development.  Although  they  are 
preliminary,  tbe  results  for  tbe  two  sets  of  conditions  are 
sufficient  for  demonstrating  the  variabihty  in  tbe  electron 
density  profile  and  its  relationship  to  corresponding  variabihty 
in  optic^  intensities. 

Nighttime  Auroral  Results 

In  the  previous  section,  we  observed  that  changes  in  source 
conditions  (changes  in  solar  activity)  lead  to  changes  in  emis¬ 
sion  brightness  without  noticeably  affecting  relative  emission 
strengths.  It  is,  thus,  the  absolute  intensities  of  optical  features 
that  concern  us  in  the  daytime  ionosphere  with  regard  to 
monitoring  variations  in  the  source  spectrum.  The  situation  is 
noticeably  different  in  the  auroral  ionosphere,  where  a  wide 
range  of  altitude  profiles  can  exist  in  volume  excitation  and 
emission  rates  due  to  variabihty  in  tbe  hardness  of  the  inci¬ 
dent  electron  spectrum.  This  leads  to  some  significant  varia¬ 
tions  in  relative  intensities,  since  changes  in  the  above  altitude 
profiles  affect  excitation  efficiencies  and  the  extent  to  which 
pure  absorption  and  quenching  operate  We  shall  now  ex¬ 
amine  such  variations  for  the  purpose  of  determining  how  well 
one  can  hope  to  characterize  the  source  spectrum  from  optical 
measurements. 

Intensities  have  been  calculated  for  viewing  in  the  nadir 
direction  from  above  the  emitting  region.  The  incident  spec¬ 
trum  is  described  by  a  Maxwelhan  distribution  having  the 
form 

♦oC  M  )  *  I I  *«p(  -  £/£m  )c/cm^-s-eV-sr 

where  £  and  ft  arc,  respectively,  the  electron  energy  and 
cosine  of  the  pitch  angle,  is  the  characteristic  energv,  and 
Qsr  is  the  enerey  flux.  Calculations  have  been  performed  for 
Qm""  1  crg/cm -s  and  for  £*,  values  over  the  range  from  0.5 
to  5  keV.  Figures  6a  and  6b  show  the  calculated  intensities  of 
the  features  under  consideration  as  functions  of  £i^.  A 
Jacchia  1977  model  atmosphere  was  used  for  all  calculations. 
We  observe  varying  amounts  of  decrea.se  for  all  features 
except  N2*  3914  and  LBH  1273  A.  The  decreases  in  Fig.  6a 
are  caused  primarily  by  the  increasing  Mount  of  O2  pure 
absorption  Some  of  the  decrease  at  1356  A  is  also  caused  by  a 
decreasing  excitation  efficiency  for  atomic^  oxygen  Pure  ab¬ 
sorption  has  little  effect  at  1273  or  3914  A,  or  at  the  wave¬ 
lengths  considered  in  Fig  6b  The  decreases  exhibited  at  2762 
and  2972  A  are  caused  primarily  by_  quenching  Some  of  the 
decrease  at  2972  A,  as  at  1356  A,  is  also  caused  bv  a 
decreasing  excilau'on  effiaency.  The  slates  from  which  these 
emissions  arise  arc  Nj(A  ’2)  and  0('S),  which  arc  coupled  in 
the  sense  that  an  important  source  of 'S  excitation  comes  from 
die  reaction 
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Nj(A)-fO-N,(X)4  o('s) 
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Fig.  6  Nadir  viewing  auronl  intcmitlM  n  funcriom  o(  tkt  incUeal 
Hcctron  vpcctnim.  .Spectra  ire  repmented  by  MaxwcSIan  dbtrlbiitiaiB 
wiiere  is  the  Maxweflian  chanetcrislic  energy. 


The  more  prominent  emission  from  the  'S  stale  is  at  5577  A, 
which  we  have  excluded  because  its  nadir-observed  intensity  is 
contaminated  by  emission  and  scattering  below  the  iono¬ 
sphere.  The  3914  A  intensity  wiU  also  be  affected  by  scattering 
and  albedo,  but  in  spite  of  this  we  include  it  in  Fig.  6  since  it 
is  one  of  the  most  prominent  and  frequently  observed  uv 
features  in  the  aurora. 


ELECTRON  DENSlTYlcm'*) 

Fig.  7  Auroral  E  region  electron  densities  for  four  different  incident 
etectron  spectra  represented  by  Maxwellian  distributions. 

We  have  obtained  profiles  similar  to  those  shown  in  Fig.  6 
for  other  model  atmospheres.  In  particular,  we  have  scaled  the 
Oj  and  0  densities  in  various  ways  to  generate  these  modeb. 
As  expected,  the  most  notable  changes  have  come  from  scaling 
the  O  density.  The  01  1356  A  intensity  varies  nearly  directly 
with  the  O  scaling  factor  because  the  main  source  is  direct 
excitation  of  0<*S)  with  a  '  lall  contribution  (-  10%)  from 
dissociative  excitation  of  O^-  Little  effect  is  observed  in  the  N2 
emissions.  Thus,  this  intensity  is  potentially  valuable  for  moni¬ 
toring  either  source  conditions  or  the  O  density,  depending  on 
how  well  these  are  known  from  independent  information.  We 
should  note  that  multiple  scattering  of  the  1356  A  photons 
does  occur,  but  has  only  a  minor  effect  for  the  given  viewing 
conditions. 

The  final  results  we  wish  to  consider  arc  electron  densities 
and,  in  particular,  the  sensitivity  of  their  altitude  profiles  to 
changes  in  the  source  spectrum.  Figure  7  illustrates  this  sensi¬ 
tivity  in  which  we  show  four  profiles  for  characteristic  energies 
of  0.5,  1,  2.5,  and  5  keV.  Since  most  of  our  attention  until 
recently  has  been  directed  to  the  auroral  E  region,  we  choose 
here  to  limit  our  discussion  to  the  behavior  of  n,  under 
chemical  equilibrium  conditions.  For  all  results,  Qn,  ~  1 
erg/cm^-s  as  before.  We  do  observe  noticeable  changes  in  the 
n,  profiles  shown,  although  not  in  maximum  values.  It  is 
worth  comparing  results  such  as  these  to  results  for  other 
types  of  source  representations.  We  have  done  so  by  earring 
out  a  scries  of  calculations  for  incident  electron  spectra 
characterized  by  narrow  Gaussian  distributions.*  We  find  the 
overall  profile  to  be  similar  where  respective  characteristic 
energies  are  chosen  to  give  the  same  altitude  of  maximum 
ionization.  For  the  more  nearly  monoenergetic  Gaussian 
source,  however,  a,  peaks  somewhat  more  sharply  and  docs 
not  penetrate  as  deeply  into  the  atmosphere. 

Much  more  can  said  about  the  results  presented  in  this 
section.  We  have  limited  the  discussion,  in  part  because  details 
are  available  elsewhere,*  and  our  emphasis  is  meant  to  be  on 
the  larger  picture  of  just  how  much  potential  optical  observa¬ 
tions  hold  for  monitoring  the  ionosphere.  We  complete  this 
paper  in  the  next  section  with  some  general  comments  at  this 
level  of  consideration. 

Discussion 

Gearly,  we  have  not  answered  the  question  of  just  how 
effectively  n,  and  the  major  neutral  densities,  especially  it(0), 
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can  be  monitored  by  optical  techniques.  This  is  a  difficult 
question  which  first  requires  statements  of  the  needed  accu¬ 
racy  in  the  deduced  quantities  and  of  instrument  capability. 
Given  this  information,  we  can  then  proceed  to  incorporate 
results  of  the  type  discussed  in  this  paper  into  a  statistical 
analysis,  out  of  which  come  deduced  density  profiles  with 
associated  standard  deviations  for  specific  optical  data  sets. 
This  should  be  done  initially  under  ideal  conditions,  which 
means  that  these  deviations  would  not  reflect  uncertainties 
introduced  by  winds  and  fields. 

What  we  are  able  to  say  at  this  time  is  that  selected  iiv 
features  do  show  sufficient  variability  in  their  intensities  to 
make  it  worthwhile  to  proceed  with  the  kind  of  analysis 
discus.sed  above.  We  shall  be  continuing  our  investigation 
along  these  lines  with  an  important  part  of  this  work  directed 
to  those  complications  which  remove  us  from  ideal  modeling 
conditions.  Among  such  complications  not  previously  men¬ 
tioned  are  the  nonanalytic  nature  of  actual  auroral  electron 
spectra  and  incomplete  spectroscopic  information  (see,  e.g.. 
Refs.  32  and  33). 
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ABSTRACT 


This  papor  presents  a  progress  report  <»n  our  efforts  toward  determining  the  feasibility 
■)f  nonitorinj  f.he  global  ionosphere  using  sateinte  observations  of  auroral  and  atrgl***  '^pt.icil 
‘‘missions.  After  giving  an  overview  of  our  approach,  we  concentrate  on  the  nighttrme  mtd !  jt.  r  t  ud<» 
lornsph^'re  for  ivhich  the  most  St raighlforward  methods  are  applicable.  There  are  several  atomic 
oxygen  emission  features  wh»ch  ire  directly  related  to  ionospheric  parameters.  Our  analysis 
indfcat‘’S  thjt  the  01  n^6  A  feitare  can  be  rno<?/tored  by  existing  photonetr'  technology.  This 
•’’eiture  provides  a  direct  measure  of  the  peak  electron  density  and  total  electron  content.  The 
01  6.100  A  feiture  can  be  used  to  monitor  thn  hmght  of  the  peak,  but  a  large  background  signal 
must  be  measured  and  jubtract<“d.  Our  analysis  indicates  that  a  state-of -lhe*irl  tilting  filter 
photometer  oan  be  used  for  this  purpose.  Howev/^r,  s  mre  no  *usfrument  cnmoarabl**  to  the  one  we 
propose  hjs  «ver  been  flown  on  \  satellito,  4  dedicated  proof -of -pr inc ip )e  experiment  should  be 
performed.  The  davs’de  and  auroral  ionospheres  have  no  comparable  f'Mt'i»*<>s  which  provid**  direct 
measures  :if  the  'Inctron  d<»nsit/.  However,  the'*<»  arp  sev»ral  emission  features  wh’rh  resoond  ^0 
cb  nges  m  the  source  j/id/or  th--*  neutry)  .>tmos/)'>erp.  We  show  synthet  ic  spectra  which  illustrate 
th’s  variability  and  briefly  di.cuss  the  usefullness  of  '>o"ie  of  these  features. 


1.  iNTROOUCTIOH 


In  this  paper,  we  address  the  problem  of  remote  optical  sensing  of  the  electron  density 
profile,  n  on  a  global  scale.  We  specifically  apply  a  nighttime  model  to  a  highly  dynamic 
situation  and  also  show  the  kinds  of  UV  spectral  variability  to  be  expected  m  the  aurora’  and 
Jaytiine  'onnspheres  for  nadir  viewing  from  a  satellite.  In  the  nighttime  ionosphere,  th»  optical 
•♦missions  of  interest  arise  from  various  chemical  processes  starting  with  the  loss  of  0*  by 
either  recomb inat ion  or  chai'ge  exchange.  Features  which  we  single  out  here  are  01  1356  A  and  01 
6300  A.  The  model  discussed  m  Section  3  shows  how  these  features  may  he  us »d  to  directly  deduce 
n  .  In  ♦he  iuroral  and  daytiine  ionospheres,  fhe  optical  "missions  of  interest  arise  from  the 
interaction  of  Secondary  and  primary  electrons  with  neutrals.  Features  which  we  single  out  here 
arp  N-  '  BH  binds  and  01  l3Sf>  A.  UV  sppctra  containing  these  ind  oi'hpr  features  vary  with  the 
suurrp  and  n<Mj*ral  atmosphere  ind  therefore  convey  information  about  those  guantitips.  This 
•of  o’-maf 'on  my  N>  j-w>d  to  Specify  n  with  .ippr'ipr  i  ate  models.  tVy  issues  h*‘re  ar®  1)  h'w  well 
',M  jrt:e  and  nputr.il  density  inf 'rtnat  i?)n  may  t>e  spprified  f'^om  satellite  observed  optical  I'missions 
and  7)  how  accurately  n^  may  deduced  once  givon  this  informatim.  JUimat‘'*ly.  these  guestinns 
mr,t  be  answerpd  with  a'dedicated  experimental  pi'ogram.  In  the  meantime,  iwe  -lefine  as  best  we 
.  ,jn  ♦hp  rn'ationship  between  n^  and  UV  *pectra  using  first  principles  models  and  jvailab)-? 
rp’pvant  datj?" 

Most  of  this  paper*  ad'Iresses  the  nfghttime  midlatitude  ionosphere.  For  orientation,  we 
list  the  soorres  of  unrertalnt/  affecting  o  as  d«*dus,ed  from  optical  'tata.  They  are  1) 

.  I’lpr.iMon,  ri  lOuntitig  stati.tics,  H  dare  <ounts,  A)  hack grtnjnd ,  inrl  S)  ^he  ipplied  •rndpl. 

T'i»*  '’/).ip|  '-'troduf  ps  uniprtam’y  throiigh  1)  an  iss'imed  shape  for  n  ('hipman  ftmct’onl^  ii  in 
a-'-umeiJ  exo-.pMer  i.;  t  mnpMr  a  L  ur  t» ,  1)  .ippli^d  f.ite  i.i*»ff*i  lynts,  jnd  as''  jmmtj  ?)*  l.’p  fhe  nly 
'in  on  the  '  '0  S  '  de  ;  ignoring  incr*»asing  *mpor  t  anre  d  u  with  '  nt  rp  a*.  i  orj  iMitudel.  T^>pse 
villous  'jn<  nr  t  a  lot  i '•‘S  w’ 1  1  tie  Msi.iissed  "i  later  mns  . 
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2.  fiACKGROUNO 


Over  the  p4St  few  vear^,  we  have  been  investigating  the  relationship  between  UV 
emissions  and  the  “lert^on  density  profile,  n  (z),  Murh  of  this  effort  has  been  directed  to  the 
continuous  aurora  (Strickland  et  a1.,  1983).  ^Recently,  we  have  turned  our  attention  to  the 
midlatitude  nighttime  ionosphere  and  the  low  to  midlatitude  daytime  ionosphere.  (Strickland  et 
al.  Our  long  term  interest  is  to  establish  the  effectiveness  of  satellite  near -nadir  UV 

data  for  monitoring  n  (z)  from  the  f  region  to  the  upper  f  region  under  quiescent  plasma 
condlMons,  A  parafi?!  effort  i'*.  under  way  by  R.R.  >1eier  and  colleagues  at  the  Naval  Research 
i.abnrjtory.  To  date,  their  ^phasis  has  lieeo  on  establishing  the  usefullness  of  UV  limb  intensity 
profiles  for  monito»''ng  the  major  neutral  constituents  (N-,  0^,  and  0),  exosphere  temperature, 
and  •excitation  sources  (solar  EUV  Spectrum  and  incident  aororjl  el^^ctron  spectrum;  see,  e.g, 

Meior  vid  Anderson,  (984  and  Meier  et  al.,  1982).  The  NRL  program  effectively  compliments  our 
own,  the  latter  being  under  Air  Force  sponsorship. 


Our  program  is  a  combination  of  1)  investigating  basic  me'hanisms,  2)  determining  what 
can  be  learned  from  '•<isting  satellite,  rocket,  and  ground  based  dita,  and  3)  directing  attention 
to  systems  aspects  '>1  the  problem  such  as  instrument  definition  and  operational  requirements . 

The  ‘nvest  iqat  lo'i  of  hisic  mechanisms  has  been  done  with  first  priciples  models  applied  to  the 
auroral  ind  diytime  ionospheres.  Here,  we  begin  with  a  source  spectrum  (energetic  incident 
auroral  electr,-»ns  or  a  solar  fiJV  and  X-ray  spectrum)  and  model  neutral  density  atmosphere.  The 
electron  density  and  numerous  UV  intensities  are  then  calculated  using  transport  and  chemistry 
algorithms  within  those  first  principles  models.  These  results  ar.»  examined  tc  determine 
patterns  of  ch  mge  with  changes  to  Input  parawters  and  are  also  compared  with  data  when 
)vai Table.  Most  of  our  efforts  have  been  directed  to  the  £  and  lower  F  regions.  We  are  now 
heq'oning  to  actively  investigate  how  well  n  can  be  determined  at  higher  altitudes  using  a 
comLioation  of  UV  data,  in  situ  measurements^at  satellite  altitudes,  and  present  knowledge  of 
F-reqion  dynamics,  A  key  to  using  the  optical  techniques  on  the  aurora)  and  daytime  ionospheres 
is  iJV  spectral  variability.  Examples  of  this  variability  will  be  given  in  Sections  7  and  9  in 
the  form  of  synthetic  spectra. 

Rugardinq  it^m  21,  Our  primary  interpst  i$  in  those  Situations  where  n  and  UV  emissions 
af'e  being  ■>  inu  1  t.ioi’ous  1  y  monitorpd.  Unfortunately,  such  situations  are  rare  anS  generally  were 
not  ^ or  /•.>rrpl  it  ing  these  observations  (rocket  pxperimpn's  are  the  exception  and  are 

nontioned  bpl^w).  Thu'.,  we  f 'n<l  oursplvos  Searching,  e,q.,  for  g»*ound  based  ionosonde  or  radar 
daf  1  rotnfilntit  wi^h  safpllite  UV  -fata  such  as  those  obtained  on  iGO-T  (Bar^h  and  Schaffner, 

19/U;  ‘uMrp  jndilees,  19/7;  Prioy  and  MPipr,  19/1),  STP  77-1  (Anil.^rson  et  al.,  1976),  P78-1 
('.‘h  I'Kratnrl  1  pt  al.,  1983;  Paresrp  et  at.,  1983),  and  53-4  iHuffn^^n  et  al.,  1980).  Such  searches 
have  hepn  doMp  with  limited  surcpss  »n  the  nighttime  Ionosphere  ('^ei^r  and  Opal,  1973;  Chandra  et 
il.,  19/S).  We  arp  not  iware  of  successful  searches  in  either  the  aurora’  or  daytime 
ionospheres . 

Tne^'e  havp  been  rocket  experiments  providing  the  opportunity  to  correlate  n  and  IJV 
data.  In  '.iirnw  cases,  o  was  obtained  by  rocket  instrumentation  and  In  others  from  tfie  ground  by 
•onosonde  or  radar.  Ex^moles  of  auroral  experiments  are  a  March  1974  experiment  reported  by  Recs 
•»t  ii.  (l9//)  ind  Sharp  et  al .  (  1979)  and  a  March  1981  experiment  largely  unreported  to  date  (see 
Swj  Jer  and  Narcisi,  1981  for'  ion  data).  We  have  compared  calculated  n  and  UV  intensities  with 
iii*4  from  both  ejrpertmpnts  and  have  found  genwraijy  good  agreement  (<  differences  over  most 
i)f  ihp  tibsprv»d  al'  itude  ranqp).  An  p<ample  of  a  daytime  experiment  is  one  in  January,  1978 
reported  by  Gentieu  et  al.  (19/9)  and  Me»er  et  al.  (1980).  In  this  case,  n  was  obtained  from 
gr  iund  based  lobosnnde  data  but  has  not  been  reported  in  the  literature.  ^he  UV  optical  data 
are  well  understood  (Meier  ef  al.,  1980)  and  we  are  now  in  the  process  of  Comparing  calculated  n^ 
va'ues  with  the  measured  ones. 

As  noted  ibove,  there  are  a  limited  number  of  opportunlNes  to  examine  n  and  UV  optical 
laM  nbtii'^ed  at  the  Same  time  and  location.  Ultimately,  an  experimental  progran^must  be 
jnderi  ikmo  using  \  '.at^llite  and  an  array  of  ground  stations  for  obtaining  these  quantities  in 
w  )#»r-  Ip  pff.'rtively  tpst  the  UV  remote  sensinq  technigup.  Until  then,  the  hpst  we  can  ilo  is  to 
cinM'iup  improving  our  under  s  t  and  ing  of  basic  ‘nechanisms  through  modeling  jnd  analysis  of 
•listing  (Ilia  even  though  they  generally  form  an  incomplete  set  'n  terms  of  sources,  opulral  and 
(h.irjpd  particle  densities,  and  optical  •wnissions. 
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Regarding  item  3),  our  attention  has  been  directed  to  Instrumentat Ion  for 

observing  the  nightglow.  The  challenge  here  U  to  achieve  the  necessary  sensitivity  for 
fnonUorlnq  e<n1ss1on  rates  at  the  sub-Rayleigh  level  and  the  necessary  wavelength  resolution  to 
provide  the  needed  information  for  subtracting  backgrounds.  Photometers  and  spectrometers  are 
both  being  assessed  for  their  effectiveness  on  the  one  hand  to  monitor  the  weak  nlghtglow 
emissions  while  on  the  other  to  monitor  the  much  brighter  emissions  In  the  auroral  and  daytime 
ionospheres . 


3.  4PPLlfO  nighttime  model 


Optical  techniques  for  monitoring  the  n  are  particularly  attractive  when  jppMed  to  the 
'‘i'^httime  mIdlUitude  ionosphere.  Here,  emission^  are  direct  signatures  of  the  height  and  ^ 

magnitude  of  n  .  These  emissions  come  from  radiative  and  ion-ion  recomb  mat  ion  (0  ♦  e  *  0*;  0 
*  O'  *  0  ♦  0*)^and  from  dissociative  recombination  (Ol  *■  e  ♦  0  •-  O").  Key  features  arising  from 
the  0*  recombination  processes  are  01  13S6  A  md  01  Ij04  A.  The  single  useful  feature  fr.jm  0., 
'■er.grib  inat  ion  Is  01  fi?')0  A  (for  nadir  viewing  rondit  ions).  We  may  relate  their  intensities  ‘ 
rnsoec t i 1  y  to  the  n  rnaximum  N  F_  and  its  location  h  F.  by 

^  *'  »>  m2  m2 


Mhere  i”!  is  tfie  column  emission  rate  and  g,  and  g,  are  functions  dependent  ■  the  assumed  shape 
of  n  and  the  chemistry  (Tinsley  and  Bittencourt,  1975  and  Chandra  et  a'.,  1175).  F inure  1)  and 
lb  sSow  the  Important  functional  dependences  within  these  equations.  Here  the  shape  of  n^  is 


Figure  la.  The  dependence  of  1356  A  Intensity 
on  the  peak  electron  density,  N^F, 
The  ordinate  is  the  Square  root  oT 
the  intensity.  It  is  a  linear 
function  of  and  has  only  a 

weak  dependence  on  h^^F^ 


Fiqure  lb.  The  dependence  of  6300  A 
intensity  on  the  height  of 
the  electron  density  peak, 
h  fy-  ordinate  is  the 
ratio  of  the  square  root  of 
the  1356  A  intensity  to  the 
6300  A  intensity 
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given  by  a  'Tv^dified  Chapman  functinn,  which  generally  character i res  n  reasonably  well  over  the 
altitude  r jnqe  containing  most  of  the  total  electron  content.  The  square  root  of  the  1356  A 
intensity  is  very  nearly  proportional  to  N  and  insensitive  to  h  F..  The  ratio  of  the  square 
root  of  the  1  356  A  intensity  to  the  6300  A^'intensity  Is  a  sensUivij  function  of  h  F.  but  nearly 
independent  of  N  F-.  Note  the  loq  like  dependence  of  »^”Ti3^/AnIg3Qo  in  FigurP  ^Ib.  This 
dependence,  especially  for  h  F.  300  km,  is  favorable  for  determining  h  F.  even  in  the  presence 
of  likely  levels  of  uncertainty  in  the  intensities.  ^ 


4.  NIGHTTIME  FLECTRON  OENSITY  ANO  EMISSION  CHARACTERISTICS 


The  midlal > lude  ionnsoh*»re  undergoes  diurnal,  seasonal,  and  solar  cycle  variations  as 
IS  I'-r  w^-lar  chv’g»’S  iss''ri|t*>d  wi»h  'nagnetic  storms.  Latitudinal  variations  are  relatively 
nitdest  \t  T\i.:dip  iiitftjdps,  hut  lenmagnet  ic  disturbance  effects  ^re  strongest  in  the  high 
'at'fijjJe  riijfiP.  The  normal  entr^mes  in  peak  density  are  5  x  10  rm”^  (usually  just  before  dawn) 
«nii  ’■>  ,m'  fusiafly  fust  ifter  sunset),  with  occasional  excursions  outside  this  range  due  to 

V’  l'^a'inp*  ir  d  1 '.t ijrh.inf  es  ind  soMr  -I'.tivity.  The  height  of  the  peak  is  almost  always  above  250  km 
(j'irrifj  *he  m  jh f  in. I  (  immonly  n  es  to  100-150  km  near  midnight  (500  km  at  the  equator). 

n.^.  ,  th**  n,  f  ,  iMd  F  rpgions  .1^  romposed  of  molecular  ions,  they  tend  to  disappear 
at  Mjh^  ’1  ipid  '  Pf  -mb -eaf  .  A  small  F  layer  Igenerally  ^*ss  than  10  cm’"^)  is 

mjintiinMii  ihrvi'.hfnjt  i he  o’ jht  hy  scattered  FUV  radiation  (Bauer,  1973;  Schunk  and  Nagy,  1980), 
The  F,  r,*gion  *s  i1o«j i 'i at .'d  by  )  kod  is  maintained  by  downward  diffusion  of  0  from  higher 
altitudes,  Th.»se  t)  '  ins  are  lost  du<»  to  ■  harge  exchange  with  molecular  species  below  the  peak. 
The  peak  .■«  i.urs  near  »he  altt^  id**  where  the  time  scale  for  vertica’  transport  equals  the  chemical 
ti-ne  ->r'‘  a  (8auer^  Because  vertical  transport  is  s  >gn  If  leant  ]y  influenced  by  neutral 

■find',  the  variability  of  the  F-  peak  is  closely  related  to  variab  lity  in  the  thermospheric 
circulation,  r.tr  example,  eguatorward  directed  winds  blow  the  pU’-ma  up  along  field  lines.  This 
produces  not  only  a  higher  altitude  for  the  F_  peak  but  also  a  greater  density  at  the  peak 
b-»ca use  hi'miral  loss  rates  d^croase  rapidly  Nith  increasing  height  due  to  decreasing  neutral 
density  's^e,  »*.g.  Fvans,  19/5;  3tr()bel  and  McElroy,  19/0).  The  reverse,  of  course,  applies  for 

OOlewjrd  <iinds. 

Of  fhp  •*missi-»n  features  which  ran  be  used  for  ior.ospher ic  monitoring,  two  important 
ones  ar*  me  r.imic  •)«y'i*»n  1 -''.'s  with  wavelengths  of  6300  A  and  1.356  A.  We  have  developed  a 
r'lmpu’^r  ,  uije  wh  i  h  »h.»  intensity  of  these  features  as  they  would  be 

iibs.»r;ed  by  a  nigh  alti’ulp  satellite.  This  code  is  based  on  the  model  described  by  Tinsley  and 
8trr‘’nrnurt  {19/5)  md  Lhjndr.j  ,‘t  al.  (19/5).  Using  the  diurnal  aud  seasonal  variations 
'll  •.DI^y^»d  in  Fi.jiim  and  fhis  .  nmpiiter  rode,  we  have  calculated  th.r  correspond ing  variation  In 
the  ^nission  fpiturps.  The  intuositins  as  a  function  of  local  time  are  shown  in  Figure  3.  Note 
*hit  while  thp  6100  A  intensity  lenerally  lies  between  10  and  100  (ayleighs,  the  1356  A  intensity 
lips  between  O.OS  and  ?  Rayleighs.  The  low  end  of  the  range  of  1356  A  intensities  is  reached 
sometime  during  fhe  night  regardless  of  season.  These  values  are  the  worst  case,  solar  minimum, 
and  will  generally  be  larger  at  solar  maximum. 


socee  MINIMUM 


•  a  <•  ?o  i«  « 
local  mean  Time 


LOCAL  MEAN  TIME 


Figure  2.  Comparison  of  thp  diurnal  variation  of  the  peak  electi'on  density  'tnd  height 

of  the  pnak  (h^F^)  uver  Millstone  Hill  for  sunmer,  winter,  and  equinox  near  solar 
minimum  (After  Evans,  1967.) 
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LOCAL  time 

Figure  3.  Local  tinie  variations  of  6300  A  and  13S6  A  emissions  for  summer,  winter,  and  equinox 
conditions.  The  calculated  intensities  are  based  on  the  electron  density  information 
in  Figure  2. 


A  representative  latitudinal  variation  of  the  6300  A  and  1356  A  intensities  is  shown  in 
Figure  4.  These  were  calculated  using  a  representative  ionospheric  variation  shown  in  Figure  5. 
(This  figure  is  a  synthesis  of  data  from  the  Alouette  and  ISIS  satellites  as  reported  by  MuMrew, 
1965,  Nelms  and  Lockwood,  196T,  and  Chan  and  Colin,  1969).  For  this  particular  example,  the  6300 
A  intensity  falls  between  1  and  10  Rayleighs  while  the  1356  A  intensity  has  bout  the  same  range 
as  in  Figure  3.  Thus  a  conservatively  designed  optical  system  intended  for  use  as  an  Ionospheric 
monitor  must  be  capable  of  detecting  as  little  as  1  Rayleigh  of  6300  and  0.05  Rayleigh  of  1356  A. 


GCOMAGNETiC  LATiruDE 

Figure  4.  Latitudinal  variallon  of  6300  and  1356  A  emissions  for  conditions  representative  of  the 
post-midnight,  winter  ionosphere  near  solar  maxiiruro.  The  calculated  intensities  are 
based  on  the  electron  densities  of  Figure  5. 
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figure  5.  latitudinal  variation  of  peak  electron  density  IN  and  the  height  of  the  peak  (h  F  0, 

The  pint  is  representative  of  ionospheric  conditions  after  midnight,  in  winter,  and'" 
"ear  solar  maximum.  (Nee  text  for  .ources.) 
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5.  factors  which  affect  the  accuracy  of  the  ocouceo 
IONOSPHERIC  PARAHETERS 


The  prohlem  of  mnn itor inq  the  Ionosphere  using  optical  emissions  may  be  separated  Into 
two  parts:  (1)  the  me isurement  of  the  intensities  of  the  optical  emissions  and  I?)  the 
ralculation  of  the  electron  density  profile  from  the  measured  intensities.  The  accuracy  of  the 
final  result  depends  on  the  uncertaint  ies  associated  with  each  step.  The  uncertainty  of  the 
first  step  is  controlled  by  the  characteristics  of  the  measuring  instrument,  the  characteristics 
of  the  emissions  being  measured,  and  the  statistical  properties  of  the  measurement  process.  The 
accuracy  of  the  second  step  depends  on  the  accuracy  of  the  model  used  to  connect  optical 
intensities  to  the  electron  density  as  well  as  on  the  uncertainties  of  the  parameters  used  in 
that  model . 

The  characteristics  of  the  instrument,  such  as  sensitivity,  affect  the  accuracy  of  the 
measurement  because  they  are  imperfectly  knoi^.  The  major  contributor  to  uncertainty  is  the 
instrument  sensitivity  which  is  determined  by  calibration  with  a  standard  light  source.  The 
leometric  factors  (aperture  and  field  of  view}  usually  contribute  very  little  uncertainty.  For 
eiample,  Huffman  et  al.  (080)  cite  uncerta  iot  ies  of  ?*,  It,  and  7-l2t  for  aperture,  field  of 
ilew,  and  sensitivity  of  the  optical  instruments  on  the  S3-4  spacecraft. 

the  m«jsuremeot  process  consists  of  counting  photons  with  an  efficiency  determined  by 
the  instrument  sensitivity.  The  number  of  counts  recorded  during  a  given  integration  time  has  a 
statistical  uncertainty  which  is  inversely  proportional  to  its  own  sguare  root.  Increasing  the 
integration  time  reduces  this  uncertainty.  In  addition,  the  instrument  has  a  dark  counting  rate 
caused  hy  energetic  particles  in  the  space  environment.  The  dark  count  must  be  subtracted  from 
the  nominal  count  and  has  its  own  statistical  uncertainty.  The  previously  clt“djS3-4  instruments 
had  dark  count  cates  of  0. 5-0.8  counts  sec*'  in  the  laboFatory,  4-7  counts  sec'  in 
spice  at  low  and  middle  latitudes,  and  25-50  counts  sec'  in  the  polar  regions. 

In  the  ultraviolet  part  of  the  spectrum,  extraneous  background  emissions  are  usually  not 
a  problem.  (See  below  for  exceptions.!  In  the  visible  part  of  the  spectrum,  however,  there  are 
iiiirmaily  two  s.iurres  of  background  .vnissions.  (Again,  see  below  for  exceptional  cases.)  The 
fic.t  t)f  these  is  due  to  extraneous  airglow  features  which  are  unrelated  to  ionospheric 
processes.  For  the  0(  hlDO  A  line,  this  means  one  or  more  lines  of  the  OH  (9-3)  band.  The  other 
source  of  ba>k'jround  is  light  reflected  from  the  ground  and  clouds  or  scattered  by  aerosols  in 
*he  lower  atmosphere.  This  includes  the  continuum  from  astronomical  sources  (moon,  stars,  and 
the  /ndiici!  lignt),  the  OH  lines,  and  'he  01  6)00  A  line  itself.  The  OH  contaminat  ion  ran  be 
moni'ru.id  (>v  jbsgrv  mq  neirhy,  isolated  OH  lines.  The  reflected  astronomical  light  can  also  be 
monitored  hy  oltserving  at  nearby  wavelengths  between  airglow  features.  Furthermore,  if  the 
inri-jent  intensity  of  (he  astronomical  source  is  known,  the  effective  albedo  can  be  calculated 
and  used  to  subtract  the  reflected  01  5300  A  intensity. 

During  magnetic  storms,  energetic  oxygen  atoms  and  ions  are  precipitated  at  midlatitudes 
from  'he  ring  ,,urcent  (Rnhrhaugh  et  il,,  1983).  This  precipitation  can  stimulate  as  much  as  20 
kR  of  01  6300  4  and  100  R  of  01  77/4  A  (equivalent  to  140  R  of  01  1356  A).  This  is  sufficient  to 
comoletely  ihscgrn  the  contribution  from  the  ambient  ionosphere.  Torr  et  al .  (1976)  reported 
ihat  t'ven  under  magnetically  guiet  conditions,  there  are  times  when  energetic  electrons  and 
protons  precipitaie  at  mid  I  at i t udes .  Using  the  electron  spectra  which  they  reported,  we  have 
cilcuiited  intensities  of  1.5  R  of  01  6300  A  and  0.04  R  of  01  1356  A.  These  are  comparable  to  the 
minimm  intensities  which  we  nvooct  to  see  from  the  nighttime  ionosphere.  Huffman  et  al .  (1880) 
jbserved  N,  '.BH  band  •■mission  at  mid lat i tudes ,  most  frequently  during  northern  summer.  The  cause 
)f  'his  emtsMon  is  'inknown  (Heier  and  Conway,  1983)  and  it  tends  to  obscure  both  the  01  1355  A 
and  the  01  1304  A  linns.  Note  that  none  of  these  events  causes  contamination  of  the  01  911  A 

lo'ij/itinn  coiitin.ium.  Although  it  is  only  half  as  bright  as  01  1  356  A,  it  could  also  be  used  as 

noMitiir  ,}f  .lU’f.tron  d#»n^ity, 

Qnrp  ibp  rpiPv.int  tnrpn^lties  (such  01  1^56  A  ^nd  01  6300  A)  have  been  determined, 
niy  ’■«  us'hI  Id  '  .ilrulitp  the  peak  electron  density  and  its  height  usmq  the  method  dpscrtbed 

»D  ’‘rt.tion  1.  w  (n  .addition  to  the  iinrerla  »nt  (es  associ-lted  with  the  mtensit>es,  there  are 

xnfert  !  int  ie»  introduted  by  the  caUuUtion  method  itselt,  first,  the  rnethod  requires  values  for 
1  Dumber  (if  pK 'meters  ivhKh  are  not  known  exactly.  These  include  enosoheric  temperature, 

'leutril  dfROsit'es,  ind  various  rates.  Second,  the  method  depends  ryn  an  >ssumed  shape 

for  the  .’leetrnn  density  profile  (the  awdlfjed  Chapman  ftmction).  linder  first  nrcumstmces,  this 
•;h<pe  will  differ  from  the  true  profile  althouqh  It  generally  reproduces  the  peak  quite  well. 

(')Oe  Section  6  for  an  example.) 
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In  order  to  assess  the  feasibility  of  oslnq  optical  emissions  to  monitor  the  nighttime 
Ionosphere*  we  have  analysed  the  capabilities  of  the  S3-4  photometer  with  its  1340  4  filter 
(Huffman  et  al.,  1980),  Because  the  full  width  at  half  maximum  of  this  filter  is  151  A,  the 
photometer  will  measure  both  01  1304  A  and  01  1156  A.  Because  the  01  1304  A  line  is  subject 
to  multiple  scattering,  there  is  greater  uncertainty  in  relating  its  intensity  to  electron 
density  than  there  is  for  01  13S6  A.  For  a  given  1356  A  intensity,  we  have  estimated  the  1304  A 
contribution  to  the  photometer  signal  using  the  multiple  scattering  analysis  of  Slricicland  and 
Anderson  (1977).  The  uncertainty  associated  with  this  estimate  has  not  been  included  in  our 
ru»'^‘»nt  analysis.  Taking  into  account  all  other  sources  of  uncertainty  (calibration,  geometric 
factors,  counting  statistics,  and  dark  count)  we  have  calculated  the  uncertainty  in  the  measured 
intensity  (1304  A  and  1356  A)  as  a  function  of  integration  time  for  several  values  of  the  01  1356 
A  intensity.  We  have  done  the  iame  thuig  for  peak  electron  density,  taking  into  account 
uncertainties  in  model  parameters.  These  results  are  shown  in  Figure  6.  We  did  not  take  into 
account  the  error  caused  by  assuming  a  Chapman  layer  shape  for  the  electron  density  profile.  It 
IS  clear  that  this  photometer  is  capable  of  observing  even  the  lowest  ‘expected  intensity  of  01 
1356  A  with  reasonable  integration  times. 


10  lo'  lo’ 
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Figure  6.  Uncertainties  in  4tiI.^--  and  M  F-  as  functions  of  integration  time  for  various 
values  of  the  01  l35b'^A”intensTty.  The  1340  A  photometer  on  the  53-4  satellite 
(Huffman  et  al.,  1980)  was  used  as  a  model  for  this  analysis. 


No  6300  A  photometer  designed  to  allow  accurate  subtraction  of  the  background  has  ever 
been  flown  on  a  spacecraft.  We  believe  that  a  state-of-the-art  tilting  filter  photometer  could 
do  the  job.  However,  our  analysis  of  such  a  hypothetiral  system  indicates  that  under  the  least 
fivorible  conditions  (full  mo*''n,  fresh  snow,  and  I  9  of  01  6100  A),  integration  times  of  60 
seconds  may  be  required.  For  a  system  located  at  800  km  altitude,  this  means  that  the  spatial 
resolution  is  i.6''of  latitude.  (Note,  however,  that  this  emission  intensity  is  reguir‘»d  only  to 
iibtain  the  height  of  the  peak.  The  peak  density  can  be  determined  on  a  much  finer  ‘cale  fmm  the 
1366  A  inten^tty  alone.)  To  i:rove  the  f*»asibility  of  this  system,  a  oecially  designed 
instrument  needs  to  be  flown  o  that  both  instrument  performance  and  the  background  problem  can 
be  properly  assessed. 
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6.  APPLICATION  OF  MODEL  tJNOER  STRONGLY  DYNAMICAL  CONDITIONS 

The  contour  plot^  of  Figure  7  Mlu^tr^te  a  major  challenge  to  any  ab-lnitlo  modeled  n 
calculation;  a  corrwon  major  nighttime  transport  event.  During  summer,  the  lower  mi d- 1  at i tudel 
show  significant  decreases  in  the  height  and  magnitude  of  n  on  many  but  not  all  nights  near 
midnight.  The  contour  plots  for  the  three  nights  of  data  shSen  in  Figure  7  are  illustrative  of 
the  normal  night-to-night  variability.  This  variability  leads  to  departures  of  n  from 
state-of -the-art  models  that  ran  be  a  factor  of  5-101  Using  the  observed  n  at  OlOO  and  0?30 
local  time  on  ?5  June,  we  have  calculated  the  intensities  of  01  1356  A  jnd  DI  6300  A,  Then, 
using  the  algorithm  of  Tinsley  and  Bitlencourt  (19/5),  we  calculated  n  corr.^sponding  to  those 
intensities.  (In  this  algorithm  the  shap^  of  ne  is  assumed  to  be  a  mo8if»ed  Chapman  function.) 
The  inferred  n  's  arw  compared  to  the  observed  n  in  Figure  8.  In  contrast  to  conventional 
ix)dels  that  haSe  been  in  error  by  about  a  factor**of  8,  our  model  (based  on  optical  emission 
lines)  agrees  within  about  .Y3I,  a  striking  improvement.  Note  that  tn/s  simulation  does  opt 
include  the  effects  of  instrumental  uncertainties  or  uncertainties  in  model  parameters. 
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ri'i'ife  7.  Contours  in  log.«  of  electron  concentrations  (electrons  cm  ) ,  on  an  altitude  versus  local 
time  geid,  showing  patterns  Of  behavior  in  the  nighttime  mid-latitude  ionosphere  over 
Arecibo  at  Ifl.J  N  geographic  (?9  N  magnetic)  latitude.  Eguatorward  neutral  winds  are 
important  in  maintaining  the  nighttime  ionosphere,  as  in  the  upper  half  of  the  figure. 
Transient  poleward  wind  components  of  unknown  origin  generally  siM?ep  across  the  entire  mid- 
latitude  sector  of  the  globe,  manifested  here  by  rapid  order-of-magnitude  drops  in  topside 
plasma  concentrations  and  transient  bottomside  increases  of  two  to  three  orders  of  magni¬ 
tude  (from  Hanson  and  Carlson,  1977). 

AACCieO  CMT*  IJIME  2S.  1969) 


Figure  8,  Calculated  and  observed  n  at  0100  .md  at  0?30  local  tune  on  ?5  June  1968.  The  observed 
n  fomes  from  Figure  7.  The  calculated  is  based  on  1)56  A  and  6300  A  intensities 
derived  from  the  observed 


y.  the  nighttime  auroral  ionosphere 


Auroral  UV  spectra  can  *»itNfblt  dramatic  differences  from  oog  to  another  because  of  1) 
c^^mposUion  change  with  altitude,  time,  and  location  and  2)  0^  pure  absorption  with  its  strong 
variation  across  the  UV  region.  An  extreme  example  of  spectral  variation  may  be  seen  in  EUV  data 
published  by  Paresce  et  al.  (198T>.  Spectra  are  shown  for  a  day  and  a  night  aurora.  The 
in  ident  electron  spectrum  is  much  softer  for  the  day  aurora  and  leads  to  preferential  excitation 
of  0  i^ompared  to  N^.  The  result  is  a  spectrum  stronger  in  0  and  0  lines  and  weaker  In  N^,  N, 

•I'ld  N  compared  to'’a  night  auroral  CUV  spectrum.  These  spectra)  'hanges  convey  ‘nformatibn  on 
the  'ncident  electron  spectrum  md  atmospheric  composition  as  already  noted  earlier.  Strickland 
“t  al .  (1983)  discuss  m  some  detail  the  connection  between  these  quantities  and  UV  spectral 
rnanges  as  well  as  the  connection  with  n  .  Here,  we  will  limit  our  discussion  to  spectral 
variability  for  nadir  viewing  from  satelt’-le  altitudes. 

Figure  9  shows  two  synthetic  spectra  between  1000  and  1800  A  for  incident  electron 
spectra  given  by  .5  and  4  keV  M.i<well'an  distributions.  Each  distcibution  is  isotropic  in  the 
incidence  plane  over  the  downward  hemisphere  and  contains  I  erg/cm'^-s.  The  original  spectra  bjve 
been  smoothed  to  give  a  resoluton  of  30  A.  The  spectra  displayed  are  partial  spectra.  Th‘>y  do 
not  include,  e.q,,  the  optical'./  thick  features  01  1304  A,  HI  I?I6  A,  and  01  1027  A.  For  the 
given  resolution,  a  number  of  features  arf»  bl»»niled.  Twenty  I.BH  bjnds  are  present  even  though 
fewer  peaks  are  exhibited  in  their  spectral  region.  The  NI  and  N[f  features  shown  arise  from 
fli^ictron  impact  on  rather  thin  on  N.  The  original  intensities  were  calculated  using  the 
transport  algor  ithm'of  StrickUnd  et  al .  (1976)  with  a  Jacchia  (  1977)  model  atmosphere  having  an 
ex  spheric  temperature  of  1000  K. 


1000  1200  1400  1600  1800 

Wavelength  (A) 

Figure  9.  Calculated  nadir  viewing  UV  spectra  for  incident  electron  spectra  given  by  .5  and  4  keV 
Maxwell ians.  The  insert  shows  the  same  spectra  on  an  f^nlarged  vertical  scale. 


Thp  <•^h1Dlte(i  differences  at  selected  wavplnngtbs  due  lu  'hjnjus  'U  the  hudn^-,',  if  'h»' 
snijr(  e  hjvw  bpgn  disruss**d  by  Strirkl.jnd  »»t  .  ll'lHl)  jml  het'u  we  will  f»e  I"  luf.  i‘i.’ 

nsA  A  fpat  ir.9  UJl  I  A  md  l8H  1  154  41  IS  part  u  u  lar  I  y  s»»nsit'vn  to  th«*  .foctron  sp**ctr(/m 

»s  a  useful  'ndirator  of  this  quantity  provided  the  0  lens’ty  is  known.  The  BH  ..paitr'jm 
1.  r  ^O  shi'wn  fi  exhibit  nearurible  rhan-^e  for  the  elf*ctron  sperira  ronsidered.  f'or  *he  4 
ise,  thi-,  herumes  b-pressed  between  I  ISO  v«d  InOO  A  due  to  absorption. 
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Uftimately  w€  wish  to  determine  how  accurately  n  In  the  E  and  lower  F  regions  can  be 
specified  ff'om  spectra  such  as  those  shown  in  Figure  9.  8e  have  the  theoretical  relationship 
between  these  quantities  for  fixed  conditions  which  include  the  neutral  atmosphere  (N  and  NO  as 
well  as  Np,  Op,  and  0)  and  duration  of  the  precipitation.  Preliminary  analysis  has  been  carried 
out  to  assign^error  bars  to  the  deduced  n  for  some  of  the  uncertain  physical  parameters.  More 
work  is  required  before  reporting  these  results. 


a.  THi  OAYTCMe  IONOSPHERE 


Noticable  changes  in  'JV  spectral  shape  were  observed  in  the  auroral  results  in  the 
P'-ev»oiJS  ;ect»on.  This  will  not  occur  in  dayglow  UV  spectra  Cover  modest  changes  In  the  solar 
;nnifh  angle  (S/A)l  since  the  altitude  profiles  of  volumes  emission  are  relatively  fixed  compared 
to  Mftse  under  «. hanging  auroral  conditions.  What  will  change  is  overall  brightness  of  dayglow  UV 
spectra  as  solar  activity  changes.  An  exception  does  occur  if  the  0  density  profile  varies 
relative  to  that  of  N  In  this  situation,  features  such  as  01  1356  A  wi 1 1  also  vary  relative  to 
s  ich  t*nturps  as  LflH  1303  A. 

P  10  shows  spectra  similar  to  those  in  Figure  9  except  for  the  dayglow.  The 

inions»ti*»s  wer«  calculated  with  the  model  of  Strickland  and  Meier  (198?)  using  a  Jacchia  (1971) 
nod**!  with  an  ojrospheric  temperature  of  1000  K.  The  solar  /enith  angle  is  60^. 

Th*'  ippli»*d  solar  f'JV  spectra  come  from  Torr  ot  al  (1979).  Below  190  A,  we  use  the  results  of 
Donnelly  md  Pope  (1973)  which  were  applied  unchanged  to  the  cases  considered. 
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Figure  10.  Similar  to  Figure  9  excopt  for  the  dayglow. 


fh<fl  Tnrr  ,i  I .  '.clar  spectra  i<f#»ntifled  l^yfin  j  -  '’06  and  71  were  used  to  generate 
fhp  spprfra  in  Figure  If),  fhe  results  ir®  intendpd'^fo  reflect  conditions  under  high  and  low 
soiir  Kf'v'fy.  The  differences  are  c lear 1 y  mnasurpab le  and  suggest  that  monitoring  of  this 
>iinrtr3i  rpqion  Will  be  useful  f  >r  determining  the  overall  energy  content  of  that  part  of  the 

spet-trum  responsible  for  laintiining  the  ionosphere.  This  region  will  also  be  useful  for 
!p'p»inin»ng  the  amount  ot  0  present  md  with  these  two  pieces  of  information,  we  Should  be  ib1e 
1  1  ppcify  n  ro^snnable  well  (s.ay  )0%)  in  the  E  and  lower  F  regions.  We  are  attempting  to 
hett.»r  gu^nt’ffv  ^his  statement  ihrough  searches  for  randi'late  optical  and  n  data  to  which  our 
p''«d  let  ions  will  be  >  fimpared.  ^ 
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